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Polarization  analysis  was  performed  on  AT&T  photonic  switching  network, 
iteration  3,  and  several  improvements  to  its  existing  design  were  discovered. 
We  discovered  that  replacing  several  linear  retarders  with  circular  retarders 
Improved  the  off  axis  performance. 

A  high  precision  single  channel  polarlmeter  was  constructed  to  supplement 
data  obtained  from  the  imaging  polarlmeter. 


Executive  Summaiy 

The  Executive  Summary  provides  a  brief  synopsis  of  the  contents  of  AFOSR  Final  Report, 
Analysis  of  Polsurizing  Optical  Systems  for  Distal  Optical  Computing  with  Symmetric  Self 
Electrooptic  Devices,  contract  #AFOSR-89-0542. 

This  research  was  motivated  by  the  need  to  develop  appropriate  metrology  to  study  the 
propagation  of  polarized  light  in  optical  interconnects  such  as  the  AT&T  mgital  optical 
computer  architectures,  and  to  adapt  polarization  aberration  theory  and  polarization  ray 
tracing,  to  understand  new  classes  of  problems  arising  in  optical  systems  tor  digital  optical 
computing. 

An  imaging  polarimeter  has  been  built  and  calibrated  at  UAH  under  this  contract.  This 
imaging  polarimeter  measures  polarization  properties  of  a  set  of  ray  paths  through  optical 
systems  and  subsystems.  The  output  is  a  Mueller  matrix  image,  a  spatial  representation  of 
the  polarization  properties  across  a  wavefront,  typically  as  a  ninction  of  field  of  view  or 
pupil  coordinafe.  The  imping  polarimeter  provides  performance  measures  of 
polarization  components.  Because  it  provideii  polarization  images  over  a  wavefront,  the 
imaging  polarimeter  presents  an  unprecedented  view  of  what  occurs  during  the  alignment 
of  optical  and  polarization  components.  With  an  imaging  polarimeter,  one  can 
simultaneously  balance  the  transmitted  polarization  between  the  top  and  bottom,  left  and 
right  sides  of  a  wavefront  as  one  rotates,  tilts,  or  otherwise  adjusts  polarization 
components..  Using  this  technique,  O.T  alignment  error  in  a  quarterwave  retarder  has  been 
resolved.  Phase  measuring  interferometry  in  real  time  had  a  similar  impact  on  optical 
alignment.  Appendix  A  contains  a  paper,  "Imaging  pblarimeters  for  optical  system 
metrology,"  which  details  the  construction  and  operation  of  the  imaging  polarimeter. 

The  most  important  study  performed  under  this  contract  with  the  imaging  polarimeter  was 
a  comparison  between  commercially  available  polarizing  beam  splitters,  specifically  to 
determine  their  suitability  for  optied  interconnects.  Large  differences  in  performance 
were  observed  between  polarizing  beam  splitters;  most  of  the  cubes  polanzing  properties 
were  surprisingly  poor.  Appendix  B,  a  copy  of  a  paper  submitted  to  Applied  Optics, 
summarizes  this  effort. 


A  polarized  spherical  wave  exiting  a  polarizing  beamsplitter  has  an  orientation  of  linear 
polarization  which  rotates  from  top  to  bottom,  a  form  of  polarization  aberration.  This 
polarization  aberration  has  been  carefully  documented,  experimentally  and  theoretically. 
The  variation  of  orientation  of  linear  polarization  results  when  the  angle  of  incidence  lies 
outside  the  normal  plane  (the  plane  containing  the  face  normals  and  interface  normal, 
usually  a  horizontal  plane  when  incident  and  exiting  beams  are  parallel  to  an  optical  table) 
Figure  1,  on  the  next  page,  shows  the  polarization  aberration  of  a  horizontally  polarized 
spherical  wave  focused  through  a  polariring  beam  splitter.  A  fraction  of  each  ray  incident 
out  of  the  normal  plane  is  reflected,  due  to  the  rotation  of  the  plane  of  incidence  rays 
within  the  spherical  wave  make  with  the  beam  splitting  interface,  causing  a  rotation.of.the. 
linear  polarization  direction. 
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Figure  1  Polarization  aberration  introduced  into  a  conver^ng  spherical  wave&ont  transmitted  through  an 
ideal  polarizing  beam  splitter.  A  uniform  horizontally  polaruxd  (within  the  plane  of  the  p^e)  spheri^  wrve  is. 
incident.  The  transmission  axis  of  the  polarudng  beam  splitter  cube  b  aligned  with  the  p  direction  a  ray  makes 
with  the  beam  splitting  interface  which  varies  with  ray  aj^e  of  mddence,  causing  apodization  and  a  spatially 
varying  orientation  of  linear  polarization  in  the  transmittra  wavefront.  The  polarization  of  the  incident 
spherical  wave  is  illustrated  in  Figure  la  and  the  transmitted  wavefront  is  shown  in  I%ure  lb. 


Figure  2  shows  how  a  polarizing  beam  sphtter  situated  between  image  planes  at  infinite 
conjugate,  similar  to  many  optical  interconnects  configurations,  introduces  a 
non-uniformity  of  intensity  in  the  image  plane  and  a  spatial  variation  of  linear  polarization 
across  and  image.  Polarization  noii-uniiormity  in  the  image  can  reduce  the  efficiency  and 
throughput  of  optical  interconnects,  giving  rise  to  extraneous  beams  at  successive 
polarization  beam  splitters.  Tlie  spatial  variation  of  linear  polarization  states  must  be 
compensated  for  to  allow  efficient  transmission  or  reflection  at  subsequent  polarizing 
beam  splitters. 

Our  analysis  of  iteration  3,  the  latest  iteration,  of  the  AT&T  Photonic  Switching  Network, 
uncovered  ways  to  improve  the  efficiency  of  polarization  optical  intercoimects.  This  study 
proceeds  from  the  supposition  that  polarizing  beam  splitters  provide  the  most  efficient  four 
)ort  access  to  the  s-SEED,  an  assumption  supported  by  research  described  in  References 
11][12][13][143.  A  polarizing  beam  splitter  cube  has  a  spatially  varying  orientation  of 
inear  polarizaticn  eigenstates.  The  goal  of  this  research  was  to  choose  polarizing  elements 
and  to  position  them  to  transform  the  polarization  state  of  the  beams  into  the  correct 
orientations  of  linear  polarization  to  match  the  polarization  eigenstate  directions  at  each 
polarizing  beam  splitter  cube.  Existing  architectures  utilize  quarter  wave  retarders  to  do 
this.  ITie  design  is  successful  for  beams  that  propagate  within  the  normal  plane  of  the 
polariring  beam  splitter.  Outside  the  normal  plane,  however,  the  beams  are  not  rotated 
into  the  correct  linear  polarization  state  to  align  with  polarization  eigenstate  directions  at 
the  polarizing  beam  splitter.  ITie  misorientation  of  the  linear  polarization  to  the 
polarization  eigenstate  directions  of  the  polarizing  beam  splitter  cubes  causes  a  fraction  of 
the  beams  to  be  routed  incorrectly.  Our  research  shews  that  at  certain  locations  in  the 
optical  interconnect,  the  linear  quarterwave  retarder  can  be  replaced  with  a  circular 
ejuarterwave  retarder  to  manipulate  the  linear  polarization  directions  of  the  beams  to  align 
with  the  polarization  eigenstate  directions  of  the  polarizing  beam  splitters.  We  have 
suggested  replacing  several  quarterwave  linear  retarders  with  circular  retarders  and  have 
suggested  that  half  wave  retarders  be  inserted  in  certain  locations  to  manipulate  the 
polarization  aberration  so  that  it  increases  the  efficiency  of  the  optical  interconnects  for 
beams  out  of  the  normal  plane. 

Two  configurations  of  a=poiarization  beam  combining  module,  one  of  the  principle 
configurations  for  polmization  based  optical  interconnects,  are  shown  in  Figures  3  and  4. 
The  existing  confi^ation,  shown  in  Figure  3,  utilizes  a  quarterwave  retarder,  QWR2,  in 
double  pass  to  convert  the  linearly  polarized  light  transmitted  by  the  polarizing  beam 
splitter  (horizontal)  into  linearly  polarized  light  which  will  be  reflected  by  the  polarizing 
beam  splitter  (vertical).  The  orientation  of  linear  polarization  after  the  second  pass 
through  QWR2  does  not  match  the  polarization  eigenstate  directions  of  the  polarizing 
beam  splitter  for  reflection,  compare  Figures  3,ILand  3,V.  The  polarizing  beam  combining 
module  siiov/n  in  Figure  4  uiili^T;S  a  circular  retarder,  CRl,  in  place  of  QWR2,  resulting  in 
a  spatially  vaiying  polarization  eigenstates  which  match  the  polarization  eigenstate 
directions  of  the  polarizing  beam  splitter  cube,  increasing  the  efficiency  of  the  polarizing 
beam  combining  module.  Circular  retarders,  rasher  than  linear  retarders,  help  keep  light 
in  the  correct  eigenstates,  and  minimize  losses. 
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Figure  2  Polarization  variation  and  intensity  non-uniformity  introduced  mto  an  image  by  an  ideal  polarizmg 
beam  splitter.  An  ideal  polarizing  beam  splitter  is  situated  in  collimated  space  between  conjugate  image  planes. 
The  intensity  non-uniformity  and  polarization  variation  of  the  image  is  due  to  the  varying  orientation  of  p 
direction  each  object  beam  encounters  at  the  beamsplitting  interface.  The  object  polarization  directions  are 
shown  in  Figure  2a  and  the  image  polarization  directions  are  shown  in  Figtire  2b. 
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Figure  3  Graphical  polarl;  rtion  analysis  of  a  polarizing  l^am  conibining  module  which  utilizes  a  polarizing 
beam  splitter  and  a  quarterwave  linear  reta>:der.  I%U’e  I  illustrates  the  eigcnpolarization  directions  of  the 

folarizing  l)eam  splitter  culte  in  transmission,  and  figure  II  shows  the  eigenpolarization  directions  in  reflection. 

igures  ill  througn  VI  illu^itrate  the  polarization  changes  before  and  after  transmission  or  reflection  from  the 
polarizinir  beam  splitter.  The  role  of  the  quarterwave  retarder  in  this  configuration  is  to  rotate  the  orientation 
of  linear.polai'ization  of  the  beams  transmitted  from  the  polari^g  beam  splitter,  depicted  in  Figure  IV,  into  the 
eigenpolarization  direaions  of  PBSl  in  reflection,  illustrated  in  Figure  II.  The  quarterwave  retarder  does  not 
rotate  the  object  beams  polarization  directions,  for  which  y^H),  into  the  correct  orientation,  causbg  a  fraction  of 
the  beams  to  be  reflected  by  PBSl  in  second  pass. 
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beam  splitter  and  a  drcnlar  retarder.  Figure  I  illustrates  the  eigenpolarization  directions  of  the  polarizing  beam 
splitter  cube  in  transmission,  and  Figure  n  shows  the  eigenpolarization  directions  in  reflection.  Figures  Ffl 
tmough  VI  illustrate  the  polarization  changes  before  and  alter  transmission  or  reflection  from  the  polarizing 
beam  splitter.  The  role  of  the  circular  retarder  in  this  configuration  is  also  to  rotate  the  orientation  of  linear 
polarization  of  the  beams  transmitted  from  the  polarizing  beam  flitter,  dwicted  in  lugure  IV,  into  the 
eigenpolarization  directions  of  PBSl  in  refl^on,  illustrated  in  Figure  II.  The  circular  retarder  rotates  all  the 
object  beams  polarization  into  the  correct  eigenpolarization  directions  of  PBSl  in  reflection,  allowing  all  the 
energy  to  be  reflected  by  PBSl  m  second  pass. 


We  will  continue  to, develop  general  design  strategies  to  boost  the  performance  of  various 
optical  interconnects  based  on  polarization  beam  combining.  From  the  polarization 
modeling  we  would  like  to  develop  further  the  formalism  to  analyze  a  broader  class  of 
polarization  based  optical  interconnect  architectures.  The  unique  nature  of  these  systems 
require  novel  analysis  methods.  From  the  analysis  studies,  a  set  of  design  rules  for 
implementii^  polarizing  optics  in  polarization  based  optical  interconnects  will  be 
developed.  Ettective  use  of  these  polarization  configurations  requires  tools  like  the 
imaging  polarimeter  for  alignment. 

A  highly  accurate  single  channel  polarimeter  has  been  constructed  to  make  high  precision 
polarization  measurements  of  thin  films  liquid  crystals,  uniaxial  and  bi-axial  crystals  and 
polarizers.  The  single  chaimel  polarimeter  will  measure  depolarization,  the  polarization 
scrambling  of  light,  diattenuation,  retardance  and  complete  Mueller  matrices  of  samples 
with  a  single  beam.  The  single  channel  polarimeter  will  also  be  used  to  calibrate  polarizing 
opti^  for  the  Imaging  polanmeter. 

Error  analysis  due  to  polarimeter  misalignment  and  retardance  error  in  the  polarizing 
elements  has  been  included  in  the  data  reduction  algorithms.  Errors  in  the  retardance  and 
orientation  of  the  retarders  in  the  Imaging  Polarimeter  produce  large  errors  in  the 
measured  Mueller  matrix.  For  example,  misorientatipn  of  1°  in  a  retarder  can  produce  a 
30%  error  in  one  of  the  Mueller  matrix  elements.  The  data  reduction  algorithms  include 
error  correction  terms  to  compensate  for  these  errors. 

More  work  is  needed  to  increase  the  accura^  of  the  imaging  polarimeter.  Presently,  the 
accuracy  is  limited  by  the  8-bit  A-D  conversion  of  the  ccd  detector.  We  plan  to  increase 
the  accuracy  of  the  imagmg  polarimeter  by  ^grading  to  a  scientific  grade  Photometries  ccd 
detector  array  with  12-bit  A-D  conversion.  Tlie  upgrade  will  increase  the  accuracy  by  a 
factor  of  ten,  allowing  the  verification  of  polarization  aberration  theories,  measurement  of 
small  polarization  enects  such  as  depolarization,  and  accurate  measurement  of  several 
interesting  and  scientifically  valuable  samples,  optical  systems  and  subsystems. 
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1.  Background  on  Digital  Optical  Computers 

The  research  described  in  this  report  involves  the  polarization  analysis  and  metrology 
of  polarization  based  optical  interconnects,  the  foundation  of  an  architectural  approach  to 
optical  computing  pursued  by  AT&T  Bell  Labs.  AT&T  optical  computing  designs  make 
use  of  2-D  arrays  of  SEED  modulators  interconnected  in  free  space.  SEEDs  and  related 
electrooptic  devices  have  been  proposed  as  potential  components  of  free  space  digital 
optical  computers  because  of  their  bistability,  high  speed,  and  low  switching  power. 

SEEDs  or  seif  electrooptic  effect  devices  exploit  the  quanmm  confined  Stark  effect  in 
quantum  well  materials  to  produce  optical  bistability.^’^  One  version,  the  symmetric 
SEED  or  s-SEED,  contains  two  SEEDs  connected  together  across  a  fixed  voltage.  Figure  1 

8.  In  this  configuration,  the  s-SEED  is  optically  bistable  with  two  states,  (1)  SEED  1  is 
highly  reflecting  and  SEED  2  is  poorly  reflecting,  and  (2)  SEED  1  is  poorly  reflecting  and 
SEED  2  is  highly  reflecting.  The  device  has  time  sequential  gain,  such  that  the  state  of  the 
device  can  be  set  with  low  power  beams,  and  read  out  subsequently  with  high  power  beams. 
Because  the  device  operates  on  the  ratio  of  reflectances,  it  is  insensitive  to  optical  power 
supply  fluctuations  if  both  beams  are  derived  from  a  single  source.  The  s-SEED  is 
especially  promising  for  digital  optical  computing  applications  because  it  has  good 
input/output  isolation  being  time  sequential.  It  does  not  require  the  critical  biasing 
required  of  most  bistable  devices. 

Spatial  light  modulators,  optical  dynamic  memory^,  all-optical  shift  registers^,  and 
optical  set-reset  latcheo  have  been  demonstrated  using  arrays  of  s-SEEDs^^.  Optical 
systems  which  compute  or  perform  interconnections  using  arrays  of  s-SEEDs  have  been 
proposed  and  are  under  constructionll'12  The  focus  of  our  studies  have  been  on 
massively  parallel  optical  processors  constructed  from  interconnected  or  cascaded  optical 
logic  modules  (COLM).  This  method  of  interconnection  is  based  on  the  use  polarizing 
elements  and  lenses.  The  polarization/lens  based  system  allows  efficient  interconnections 
between  devices.  These  COLMs  can  be  interconnected  to  form  gates,  adders, 
programmable  logic  arrays,  and  other  components  of  conventional  electronic  computers. 
Figure  2  is  a  schematic  of  a  single  COLM.  The  section  of  the  COLM  shown  separately  in 
Figure  3,  is  used  to  clear  the  s-SEED  array.  The  COLM  reads  out  the  s-SEED  array, 
operates  on  the  information,  and  passes  data  to  another  COLM,  where  another  s-SEED 
array  is  set.  This  system  has  been  constructed  at  AT&T  Bell  Laboratories.  Operational 
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detaiis  of  this  optical  system  are  available  in  References  13,14. 

The  Cascadable  Optical  Logic  Module  depends  on  the  control  of  polarized  light  for 
its  operation.  If  careful  control  of  polarization  is  not  maintained,  separate  beams  will  not 
be  properly„combined.and.the.computer.becomes  unreliable..  The  optical  interconnects 
have  large  numbers  of  retarders  and  polarizing  beamsplitters  to  permit  separate  beams 
carrying  digital  information  to  be  combined  and  address  the  appropriate  quantum  wells  on 
thes-SEED. 

Polarization  based  optical  interconnect  architectures  suffer  from  substantial  amounts 
of  instrumental  polarization,  the  polarization  properties  associated  with  ray  paths  through 
the  system.  The  polarizing  beams  splitters  ar,  responsible  for  the  majority  of  the 
instrumental  polarization.  This  research  includes  studies  done  on  polarizing  beam  splitter 
cube  properties  to  understand  how  they  contribute  to  instrumental  polarization.  Two 
properties  of  polarizing  beam  splitters  contribute  significant  instrumental  polarization:  (1) 
the  polarizing  beam  sphtter  cube  beamsplitting  interface  often  displays  strong  angular 
dependence  of  s  and  p  transmittance  and  reflectance,  (2)  the  angular  orientation,  relative 
to  a  global  coordinate-system,  of  the  s  and  p  directions  (linear  polarization  directions) 
depends  on  ray  angle  of  incidence.  The  first  property,  the  angular  dependence  of  s  and  p 
transmittance  and  reflectance,  can  be  optimized  with  multilayer  design  techniques.  15, 16 
The,  second  property  is  due  to  the.  geometry  of  the  cube.  The  polarizing  axes  of  the 
polarizing  beam  splitter  are  aligned  with  the  p  direction  a  ray  makes  with  the  beamsplitting 
interface  in  transmission,  andlthe  s  direction  in  reflection.  The  angular  orientations  of  s 
and-.p,  however,  depend  on  the  angle  of  incidence  of  the  ray,  causing  a  angular 
nonuniformity  of  linear  polarization  direction.  Analysis  and  measurements  made  with  the 
Imaging  Polarimeter  on  polarmng  beam  splitters  are  detailed  in  "Angular  dependence  of 
Polarizihg  beam  splitters,"  submitted  to  Applied  Optics,  attached  to  the  end  of  this  report. 

Quarterwave  retarders  were  analyzed  and  measured  with  the  Imaging  properties. 
Quarti',  zero  order  quarter  wave  retarders  display  a  much  larger  angular  bandwidth  than 
the  polarizing  beani  splitters.  However,  several  properties  can  degrade  their  performance, 
such  as  crystal  axis  misorientation,  thickness  variations,  and  misalignment  of  the  fast  axis  in 
the  optical  system. 

The  first  half  of  this  repoft  analyzes  optical  interconnect  architectures,  and 
suggestions  for  improvements  to  existing  architectures.  The  second  part  describes  the 
instruraentation  developed  for  metrology  of  polariring  optical  systems,  and  measurements. 


Sections  2.1  and  2.2  provides  polarization  analysis  on  the  polarizing  components  of 
polarization  based  optical  interconnects.  Section  2.3  models  the  polarizing  properties  of 
Iteration  3  Photonic  Switching  Network,  and  makes  several  suggestions  to  improve  its  off 
axis  efficiency. 

The  major  thrust  of  this  research  has  been  in  the  development  of  polarimetric 
instrumentation  to  measure  the  polarizing  properties  of  polarization  based  imaging 
systems.  The  motivation  for  the  development  of  this  instrumentation  was  the  need  for  a 
metrology  tool  to  study  the  large  instrumental  polarization  of  polarization  based  optical 
interconnects.  Before  this  research,  the  needed  instrumentation  did  not  exist.  Section  3 
describes  the  Mueller  Matrix  Imaging  Polarimeter,  built  and  calibrated  under  this  contract, 
which  measures  polarization  for  a  set  of  ray  paths  through  an  optical  system,  i.e,,  it 
measures  the  Mueller  matrix  of  the  system  as  a  function  of  object  or  pupil  coordinate.  The 
Mueller  matrix  provides  information  about  the  diattenuation  and  retardance  over  the  pupil 
or  field  of  view. 

Section  4  introduces  a  closely  related  instrument,  also  developed  under  this  contract, 
a  Stokes  Imaging  Polarimeter  which  measures  polarization  properties  across  a  wavefront. 
The  output  of  this  instrument  supplements  information  about  a  wavefront  obtainable  by 
from  interferometry  (wavefront  error),  it  measures  the  polarization  content  of  a  wavefront. 
Both  interments  are  configured  to  operate  at  633nm  and  850nm,  the  operating  wavelength 
ofs-SEED’s. 

A  third  instrument  was  developed  to  support  the  first  two,  a  single  channel 
polarimeter.  The  function  of  this  polarimeter  is  to  align  the  polarizing  optics  of  the 
Mueller  matrix  imaging  polarimeter  and  the  Stokes  imaging  polarimeter.  Section  5 
provides  a  brief  description  of  single  channel  polarimeter  and  describes  the  alignment 
procedure  for  the  polarizing  optics. 

Section  6  discusses  the  calibration  issues  involved  in  building  the  imaging 
polarimeters.  Finally  concluding  remarks  are  provided  in  Section  7. 
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2.  Polarization  Modeling  of  Polarization  Based  Optical  Interconnects 

This  section  develops  the  field  of  view  polarization  properties  of  a  single  module  of 
the  AT«&T  Iteration  3  Photonic. Switching  Network.  First  we  describe  the  angular 
polarization  behavior  of  the  polarizing  elements,  the  polarizing  beam  splitters  and  the 
quarterwave  retarders.  Next  we  formulate  a  polarization  model  of  a  single  stage  of  the 
Photonic  Switching  Network,  and  make  suggestions  to  improve  its  efficiency.  The 
analytical  approach  used  here  can  be  applied  to  other  AT&T  polarization  based  optical 
interconnects. 

2.1  Field  of  view  effects  of  the  polarizing  beam  splitters 

This  sectioh  describes  the  field  of  view  properties  of  polari?ing  beam  splitter  cubes. 
The  work  presented  here  helps  develop  the  polarization  analysis  of  the  Photonic  Switching 
Network  in  Section  2.3.  A  more  thorough  treatment  of  polarizing  beam  splitters,  along 
with  measurements  made  with  the  Imaging  Polarimeter  are  given  in  "Angular  dependence 
of  polarizing  beam  splitter  cubes,"  submitted  to  Applied  Optics.,  in  Appendix  B. 

Ai  ideal  polarizing  beam  splitter  cube  divides  the  incident  light  into  s  and'p 
polarization  states,  transmitting  the  p  polarized  component  and  reflecting  the  s  polarized 
component.  Polarizing  beam  splitters  have  a  limited  range  of  angle  of  incidence  and 
wavelength  for  which  they  are  effective.  Polarizing  beam  slitters  sold  through  catalogs  are 
optimized  to  provide  a  large  spectral  range  within  a  small  angular  field  of  view.  Optical 
computing  architectures  typically  require  the  opposite,  i.e.,  a  wide  field  of  view  with 
quasimonochromatic  light. 

The  angular  performance  of  the  polarizing  beam  splitter  depends  on  two  factors: 

1.  The  transmittance  and  reflectance  of  s  and  p  polarization  states, 

2.  The  dependence  of  s  and  p  angular  orientations  on  angle  of  incidence. 

Research  has  been  done  to  increase  the  angular  range  for  which  p  transmittance  and 
s  reflectance  remain  high.  15,16  AT&T  Laboratories  has  designed  a  polarizing  beam 
splitter  to  provide  a  ±7°  FOV.  Tests  performed  at  AT&T  and  at  UAH  agree  with 
predictions  of  the  design.  UAH  measurements  on  the  AT&T  polarizing  beam  splitter, 
along  with  measurements  on  three  other  commercially  available  polarizing  beam  splitters 
are  presented  in  the  above  referenced  paper. 

The  geometry  of  the  polarizing  beam  splitter  causes  a  spatially  varying  orientation  of 
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linear  polarization  in  transmitted  and  reflected  wavefronts.  The  orientation  of  linear 
polarization  transmitted  from  a  polarizing  beam  splitter  dependent  on  the  ray  angle  of 
incidence.  For  a  spherical  wave  of  unpolarized  light  is  incident  on  a  polarizing  beam 
splitter,  the  reflected  and  trajismitted  beams  have  a  spatially  varying.linear  polarization. 
Since  the  orientation  of  the  pigne  of  incidence  a  ray  makes  with  the  beamsplitting  interface 
depends  on  the  angle  at  which  the  lay  is  incident  i.e.,  the  s  and  p  directions  depend  on  ray 
angle  of  incidence.  Since  the  orientation  of  linear  polarization  of  a  polarizing  beam  splitter 
aligns  with  the  p  direction  in  transmission  and  the  s  direction  in  reflection,  the  orientation 
of  linear  polarization  of  a  polarizing  beam  splitter  is  a  function  of  ray  angle  of  incidence. 
Figure  4  introduces  a  coordinate  system  for  describing  polarizing  beam  splitter  cubes.  The 
origin  is  centered  on  the  enn-jvnce  face  with  the  z  axis  normal  to  the  face.  The  beam 
splitting  interface  normal  lie;  in  the  x-z  plane  and  is  45°  from  the  z  axis.  The  direction  of 
an  incident  ray  is  specified  by  its  direction  cosines  {l,m,n)  where  I  =  cos  a,m  =  cos  3  and 
n  =  cos yand  +  m^  +  n.^  =  \ .  a ,  3 ,  and  y  are  illustrated  in  Figure  4.  For  incident  rays 
with  3  =  90° ,  the  p  direction  is  parallel  to  the  x-z  plane  (we  will  refer  to  this  as  the 
horizontal  plane)  and  the  s  direction  cligns  with  the  y-z  plane  (vertical).  For  angles  of 
incidence  in  which  3  >  90°,  the  p  plane  .'•otates  counterclockwise  with  the  respect  to  the 
horizontal,  while  for  3  <  90°  it  rotates  clockwise.  Figure  5  illustrates.the  orientation  of  s 
and  p  directions  as  functions  of  the  direction  cosines  of  the  incident  rays.  In  Figure  5,  a 
increases  from  90°  in  the  positive  x  direction  and  3  increast^s  from  90°  in  the  positive  y 
direction.  The  thin  cross  hairs  are  parallel  with  the  horizontal  and  vertical  planes  (x-z  and 
y-z  planes  in  Figure  4),  The  large  circles  indicate  the  zones  where  the  angles  of  incidence 
are  Y  =  5°  and  y  =  10°,  The  angle  (j)  is  the  angle  between  the  p  direction  and  the 
horizontal  x-z  plane. 

In  Section  2.3  we  explore  how  the  variation  of  orientation  of  linear  polarization 
degrades  the  performance  of  an  optical  interconnect.  Methods  to  correct  the  degradation 
of  the  Photonic  Switching  Network  off  axis  efficienc_y  caused  by  this  effect  are- described. 
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Figure  5  Orientation  of  the  s  and  p  directions  on  the  beamsplitting  interface  of  a 
^lanzmg  beam  splitter  plotted  as  a  function  of  beam  angle  of  incidence  on  the  cube  face. 
Die  orientation  of  the  beam  splitting  interface  is  shown  in  Figure  1.  Figures  1  and  2  share 
the  same  coordinate  system,  a  increases  from  90“  in  the  positive  x  direction.  13  increases 
from  90  in  the  positive  y,  and  the  large  circles  indicate  the  zones  where  the  angles  of 
incidence  are  7  =  5°  and  y=  10”.  Incident  rays  for  which  |3  >  90°  involve  a 
countercloc^ise  rotation  of  s-p  coordinates  and  rays  in  which  |3  <  90°  the  rotation  is 
clock\yise.  The  polarizing  axis  is  parallel  to  the  p  direction  in  transmission  and  the  s 
direction  in  reflection.  Polarizing  beam  splitters  have  a  spatially  vaiying  orientation  of 

in#»5ir  nninri'70fir\n  *  ^  j  ^ 
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2.1.1  Throughput  of  a  COLM  as  a  function  of  transmittance  and  reflectance  from  the 
polarizing  beam  splitters 

Figure  6  is  a  plot  of  the  fraction  of  light  transmitted  through  a  Cascadable  Optical 
Logic  Module.as.a  function  of  the.transmittance  of  the  polarizing  beam  splitters 
(reflectance  of  the  polarizing  beam  splitters  are  assumed  to  be  1).  For  this  analysis,  the 
transmittance  of  each  of  the  polarizing  beam  splitters  varies  by  the  same  amount.  The  light 
is  transmitted  through  polarizing  beam  splitters  4  times  and  reflects  6  times  (see  Figure  2). 
The  efficiency  of  the  COLM  strongly  depends  on  the  transmittance  of  the  polarizing  beam 
splitters.  If  the  transmittance  of  the  polarizing  beam  splitters  is  reduced  to  75%,  the 
efficiency  of  the  COLM  is  reduced  by  a  half.  The  dependence  is  nearly  exponential. 

Figure  7  shows  the  fraction  of  light  transmitted  through  a  Cascadable  Optical  Logic 
Module  as  a  function  of  the  reflectance  of  the  polarizing  beam  splitter  cubes  while  the 
transmittances  are  held  fixed  at  1.  The  COLM  throughput  is  more  sensitive  to  variation  of 
reflectance  of  the  polarizing  beam  splitters  than  transmittance  because  more  reflections 
occur.  If  the  reflectance  of  the  polarizing  beam  splitters  is  reduced  to  75%,  the  efficiency 
of  the  COLM  approaches  0. 

The  energy  lost  takes  the  wrong  direction  at  the  polarizing  beam  splitters,  the  beams 
are  transmitted  when  they  should  reflect  or  the  beams  are  reflected  when  they  should  be 
transmitted. 

2.2  Field  of  view  effects  of  the  retarders 

Several  characteristics  of  birefringent  retarders  introduce  polarization  aberration  in 
polarization  critical  systems.  First,  the  thickness  of  the  birefiringent  material  may  be  wrong 
or  wedged,  yielding  an  incorrect  or  spatially  varying  retardance.  The  crystal  axis  can  be 
misoriented  relative  to  the  faces  or  the  mount.  The  retarder  can  be  misaligned  in  the 
optical  system,  either  tilted  or  with  its  fast  axis  at  the  wrong  orientation. 

Three  other  problems  with  retarders  are  more  fundamental,  angular  field 
dependence,  diattenuation,  and  birefraction. 

The  angular  dependence  of  a  uniaxial  crystal  is  illustrated  in  Figure  8.  For  a  crystal 
axis  in  this  orientation,  two  effects  influence  its  retardance.  First,  the  difference  between 
the  extraordinary  and  ordinary  refractive  index,  or  birefringence,  decreases  quadratically 
with  increasing  field  angle.  This  property  decreases  the  amount  of  retardance  the  crystal 
displays.  The  second  effect  increases  tends  to  offset  the  first  effect  somewhat.  This  effect 
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Figure  6  Fraction  of  transmitted  light  through  a  Cascaddable  Optical  Logic  Module, 
plotted  as  a  function  of  the  transmission  coefficients  of  the  polarizing  beam  splitter  cubes. 
The  transmission  coefficients  of  the  each  of  the  polarizing  beam  splitters  were  varied  by 
the  same  amount,  while  the  reflection  coefficients  were  held  at  1. 


Figure  7  Fraction  of  transmitted  light  through  a  Cascaddable  Optical  Logic  Module, 
plotted  as  a  function  of  the  reflection  coeffiaents  of  the  polarizing  beam  splitter  cubes. 
The  reflection  coefficients  of  the  each  of  the  polarizing  beam  splitters  were  varied  by  the 
same  amount,  while  the  transmission  coefficients  were  held  at  1. 
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Variation  of  Retardance 
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Figure  8  Variation  of  retardance  in  a  uniaxial  crystal.  The  birefringence  varies  with  angle 
of  incidence  according  to  the  index  ellipsoid.  The  retardance  varies  because  of  the  change 
in  the  birefringence  and  the  varying  path  length  with  angle  of  incidence.  The  magnitude 
and  orientation  of  linear  retardance  varies  quadratically  with  angle  of  incidence.*" 


increases  the  retardance  simply  due  to  the  increase  in  the  physical  path  length  with  angle. 

Figure  9  shows  the  variation  in  retardance  with  field  angle  in  quartz  linear  retarder. 
The  crystal  axis  is  aligned  to  the  y  axis.  Shown  is  one  quadrant  with  the  lower  left 
corresponding  to  a  ray  at  normal  incidence,  and  the  angle  of  incidence  increasing  along  the 
.V  and  y  axes.  By  symmetry,  it  is  only  necessary  to  view  the  angular  dependence  in  one 
quadrant.  For  an  ideal  quarter  wave  retarder,  about  a  1%  variation  of  retardance  is  found 
over  a  10  degree  field  of  view.  The  field  of  view  dependence  of  quartz  quarter  wave 
retarders  should  not  contribute  significant  errors  into  AT&T  digital  optical  computing 
architectures. 

Figure  10  contains  a  plot  of  the  Mueller  matrix  as  a  function  of  angle  of  incidence, 
taken  with  the  UAH  Imaging  Polarimeter  over  a  ±6“  FO  V,  of  a  visible  zero-order  quartz 
quarter  wave  retarder  at  633nm.  Each  of  the  Mueller  matrix  elements  are  shown  and 
labeled  except  m  1 1,  since  we  normalize  by  m ,,  (see  Section  3.2).  The  data  of  interest  lie  in 
the  central  circular  portions.  The  center  of  each  of  the  plots  represents  an  axial  ray,  and 
the  angle  of  incidence  increases  in  the  radial  direction  out  to  6°  in  each  direction.  The 
retarder  displays  a  uniform  retardance,  within  1%,  over  the  6  “field  of  view. 

2.2.1  Orientation  error  in  the  quarter  wave  retarders 

An  analytical  study  has  been  done  to  determine  the  loss  of  efficiency  of  the  COLM 
due  to  misorientation  of  the  quarter  wave  retarders.  Figure  11  shows  the  percent 
transmission  of  light  that  remains  in  the  correct  polarization  state  as  a  function  of 
orientation  error  of  the  quarter  wave  retarders.  Each  retarder  received  the  same  amount 
of  orientation  error.  The  percent  transmission  refers  to  the  percentage  of  light  emitted 
from  the  laser  diodes  that  travels  through  the  entire  COLM  to  the  next  COLM,  assuming 
ideal  quarterwave  retarders  and  polarizing  beam  splitter  cubes. 

2.3  Polarization  analysis  of  Iteration  3  Photonic  Switching  Network 

Iteration  3  Photonic  switching  Network  utilizes  a  polarization  based  imaging  system 
to  interconnect  2-D  arrays  of  optical  logic  devices.  It  uses  polarizing  beam  splitters, 
quarter  wave  retarders  and  patterned  reflectors  to  interlace  input  beams  onto  a  modulator 
and  separate  output  beams  reflected  from  the  modulators.  Each  beam  must  pass  through 
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Figure  1 1  Fraction  of  transmitted  intensity  through  a  Cascadable  Optical  Logic  Module 
plotted  as  a  function  of  orientation  error  of  the  quarter  wave  retarders.  Each  oi  the 
quarter  wave  retarders  recieved  the  same  amount  of  orientation  error. 
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polarizing  beam  splitters  6  times,  quarter  wave  retarders  10  times  and  change  polarization 
state  14  times. 

Figure  12  shows  the  configuration  of  the  Photonic  Switching  Network.  The  line 
indicates  the  path  of  light  propagation.  The  desired  polarization  changes  of  the  beam  are 
indicated  on  the  drawing.  H  indicates  horizontal  polarization  (within  the  plane  of  the 
page),  V  indicates  vertical  polarization,  and  R  and  L  indicate  right  and  left  circular 
polarization.  The  path  analyzed  is  as  follows:  from  the  binary  phase  grating  the  beam 
passes  through  a  space  of  the  patterned  reflector,  through  the  polarizing  beam 
splitter/quarter  wave  retarder  combination  onto  the  s-SEED  in  module  1.  After  reflecting 
off  the  s-SEED  it  reflects  off  the  polarizing  beam  splitter  toward  the  interconnection  unit. 
The  beam  reflects  off  the  plane  mirror  into  the  interconnection  unit  and  the  polarizing 
beam  splitter  divides  the  beam  into  two  equal  amplitudes  in  orthogonal  directions.  Half 
the  amplitude  passes  through  the  polarizing  beam  splitter  to  the  plane  mirror  and  half 
reflects  toward  the  retro-reflector.  The  beam  which  reflects  to  the  retro-reflector  is  flipped 
about  the  vertex  of  the  retro-reflector  and  passes  back  through  the  polarizing  beam  splitter 
into  the  next  module.  Both  beams  either  transmit  or  reflect  from  the  polarizing 
beamsplitter  and  reflect  from  the  patterned  mirrors  shown  and  are  imaged  onto  the 
s-SEED  in  module  2. 

Our  analysis  was  done  over  a  ±10°  field  of  view.  The  assumptions  are  as  follows: 

1.  The  beams  are  uniformly  circularly  polarized  at  the  first  patterned  reflector. 

2.  The  quarter  wave  retarders  have  exactly  a  quarter  wave  of  retardance. 

3.  The  retardance  varies  insignificantly  over  the  ±10[EQUATION]  field  of  view.  (Section 
2.2  supports  this  assumption). 

4.  The  polarizing  beam  splitters  transmit  100%  p  polarization  and  reflect  100%  s 
polarization  for  all  angles  of  incidence.  This  assumption  is  also  made  to  simplify 
interpretation  of  results  and  can  be  readily  generalized  to  the  non-ideal  case. 

Figure  13  shows  a  chief  ray  path  through  this  portion  of  the  Photonic  switching 
Network  (Iteration  3),  The  patterned  reflectors,  mirrors,  retro-reflectors  and  modulators 
are  situated  in  image  planes.  The  polarizing  beam  splitters  and  quarter  wave  retarders  are 
positioned  at  or  near  the  pupils.  Since  the  system  is  doubly  telecentric  and  the  image 
planes  are  at  infinite  conjugates,  the  beams  are  collimated  when  they  reach  the  polarizing 
beam  splitters  and  quarter  wave  retarders.  It  is  therefore  sufficient  to  trace  a  set  of  chief 
rays  through  the  system  to  characterize  the  polarizing  behavior  associated  with  each  beam. 
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The  polarization  ray  trace  begins  at  the  first  pattern  reflector  and  ends  at  the  s-SEED 
array  in  Module  2.  The  results  are  all  given  in  terms  of  the  object  coordinates  d^^fined  at  the 
first  patterned  reflector. 

The  polarization  changes  that  take  place  in  the  Photonic  Switching  Network  are 
formulated  in  terms  of  the  Jones  Calculus,  and  are  illustrated  via  a  series  of  graphical  maps 
of  polarization.  The  graphical  plots  of  polarization  depict  the  polarization  state  of  the 
chief  rays.  This  approach  provides  an  intuitive  picture  of  the  changes  of  polarization  of  off 
axis  rays. 

Figure  14  illustrates  the  path  from  the  first  patterned  reflector  PRl,  through  the  first 
polarizing  beam  splitter  PBSl  to  the  s-SEED  which  reflects  the  light  back  to  PBSl  which  in 
turn  reflects  it  toward  the  Interconnection  Unit.  The  coordinate  system  for  the  Chief  rays 
are  defined  at  PRl:  the  positive  x  direction  is  shown  at  PRl  and  y  points  out  of  che  page. 
Like  positions  on  each  of  the  plots  I- VI,  correspond  to  the  same  chief  ray.  The  positive  x 
direction  on  plots  I- VI  corresponds  to  the  Chief  ray  path  shown  on  the  drawing.  The  linear 
polarization  states  of  the  Chief  rays  at  the  polarizing  beam  splitter  interface  are  shown 
because  it  is  easier  to  interpret  linear  polarization  states  than  circular  and  elliptical  states 
found  at  the  image  planes. 

The  circularly  polarized  beams  at  PRl  are  converted  to  uniform  horizontal  linear 
polarization  by  QWRl  just  before  PBSl  (Figure  14-III).  Figure  14-1  illustrates  the 
variation  of  orientation  of  p  directions  the  chief  rays  make  with  the  beamsplitting  interface 
of  PBSl,  shown  as  a  function  object  coordinate  of  the  Chief  rays.  Chief  rays  from  object 
points  for  which  y  >  0  encounter  p  directions  rotated  in  the  clockwise  direction  and  when  y 
<  0  the  rotation  is  counter-clockwise.  The  amount  of  rotation  is  denoted  0(  x ,  y  ^  and  is  a 
function  of  object  coordinates.  For  an  ideal  polarizing  beam  splitter,  only  light  polarized  in 
the  p  direction  is  transmitted.  To  determine  the  fraction  of  the  incident  horizontally 
polarized  beams  which  are  transmitted  by  PBSl,  the  incident  beams  polarization  state  must 
be  decomposed  into  the  local  s  and  p  directions  at  the  PBSl  beam  splitting  interface.  For 
Chief  rays  in  which  y  =  0,  the  horizontal  direction  aligns  with  the  p  direction,  and  100%  of 
these  beams  are  transmitted.  For  Chief  rays  in  which  y  5^  0,  a  fraction  of  the  horizontal 
polarization  lies  in  the  s  plane  and  is  reflected.  Because  a  fraction  of  intensity  is  reflected 
out  of  the  incident  horizontally  polarized  beams,  the  transmitted  beams  linear  polarization 
direction  is  rotated  by  an  amount  0(.v ,  y )  to  align  with  the  p  directions  of  the  PBSl 
interface  (compare  Figures  14-1  and  14-IV). 
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III.  Polarization  IV.  Polarization  V.  Polarization  VI.  Polarization 
before  PBSl  after  PBSl  before  PBSl  after  PBSl 

first  pass  first  pass  second  pass  Second  pass 


QWR2  converts  the  linearly  polarized  beams  into  right  circular  and  elliptically 
polarized  beams.  After  reflection,  the  beams  are  left  circular  and  elliptically  polarized. 

The  second  pass  through  QWR2  converts  the  left  circular  and  elliptical  beams  back  into 
linear  states.  In  double  p^s,  the  quarter  wave  retarder  acts  like  a  half  wave  retarder. 
Figure  15  shows  how  a  quarter  wave  retarder  rotates  linear  polarization  in  double  pass. 

The  linear  polarization  of  the  entering  and  exiting  beams  are  not  orthogonal  when  ys^O,  but 
have  an  angular  separation  of  2s.  Because  the  s  and  p  directions  are  orthogonal  on  a 
polarizing  beam  splitter  cube,  the  entering  and  exiting  beams  polarization  directions  must 
be  orthogonal  for  optimal  efficiency.  Because  the  linear  polarization  directions  before  the 
first  pass  through  QWR2  and  after  the  second  pass  through  QWR2  are  not  orthogonal,  i.e., 
separated  by  2s,  the  polarizing  beam  combining  module  does  not  operate  optimally  with  a 
quarter  wave  retarder.  Figure  16  shows  the  fraction  of  horizontally  polarized  light  which 
passes  through  the  polarization  beam  combining  module  over  a  ±10°  FOV. 

Figure  17  shows  the  action  of  a  quarter  wave  circular  retarder  on  incident  linearly 
polarized  light.  In  double  pass,  the  quarter  wave  circular  retarder  is  equivalent  to  a  half 
wave  circular  retarder.  The  half  wave  circular  retarder  rotates  incident  linear  polarization 
states  90°,  independent  of  incident  linear  polarization  state  orientation.  The  use  of  a 
circular  retarder  in  a  polarization  beam  combining  scheme  allows  a  more  effective  means 
of  orienting  the  linear  polarization  state  into  the  correct  direction  for  the  second  pass  (the 
reflection)  off  the  polarizing  beam  splitter  cube.  Figure  18  shows  the  polarization  states 
before  and  after  each  pass  through  the  polarizing  beam  splitter.  Note  that  the  orientations 
of  linear  polarization  for  each  pass  through  PBSl  align  with  the  s  and  p  directions  of  PBSl. 
Figure  19  shows  the  fraction  of  horizontally  polarized  light  passes  through  the  polarization 
beam  combining  module  over  a  ±10°  FOV. 

Figure  20  gives  the  difference  in  the  fraction  of  light  that  passes  through  the 
polarization  beam  combing  module  between  using  a  quarter  wave  linear  retarder  and  a 
quarter  wave  circular  retarder.  The  Jones  calculus  of  the  above  situation  is  provided  in 
Figures  21  and  22. 

The  configuration  in  this  situation  occurs  two  other  times  in  the  path  of  our  analysis. 
It  occurs  in  the  plane  mirror  path  of  the  ICU  and  it  also  occurs  module  2  (refer  to  Figure 
13  ).  Figure  23  shows  the  locations  where  it  would  be  beneficial  to  replace  the  quarter 
wave  linear  retarders  with  quarter  wave  circular  retarders. 

The  polarization  ray  trace  proceeds  in  Figure  24.  This  Figure  represents  the  path 
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Action  of  a  Quarter  Wave  Retarder  on 
Linear  polarization  in  double  pass 
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Figure  15 


Emciency  of  a  Polarization  Beam 
Combining  Module  for  incident 
horizontally  polarized  light 

(utilizing  quarter  wave  retarders) 


*  Plotted  as  a  function  of  object  coordinate,  defined  at  PRl 


Action  of  a  Circular  Retarder  on  Linear 
polarization  in  double  pass 


Quarter  Wave  Circular  Retarder  in  double  pass 
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Figim  17 


Efficiency  of  a  Polarization  Beam 
Combining  Module  for  incident 
horizontally  polarized  light 


(utilizing  circular  retarders) 


Efficiency  difference  between  a 
Polarizing  Beam  Combining  Module 
which  utilizes  circular  retarders  and 
one  which  utilizes  quarter  wave 

retarders 


*  Plotted  as  a  function  of  object  coordinate,  defined  at  PRl 
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Figure  20 


Jones  calculus  for  a  Polarization 
Beam  Combining  Module 

(utilizing  quarter  wave  retarders) 
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Figure  21 


Jones  calculus  for  a  Polarization 
Beam  Combining  Module 

(utilizing  circular  retarders) 
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Figure  22 
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III.  Polarization  IV.  Polarization  V.  Polarization  VI.  Polarization 
before  PBSl  after  PBSl  before  PBS2  after  PBS2 


from  the  s-SEED  in  module  1  to  the  ICU  with  existing  polarizing  optics.  QWR3  converts 
the  linear  polarization  states  reflected  from  PBSl  to  circular  and  elliptical  states.  The 
orientation  of  their  major  axes  are  flipped  by  the  plane  mirror  before  the  ICU  as  shown. 

The  orientations  of  the  polarization  ellipses  are  such  that  more  power  is  transmitted  by 
PBS2  than  reflected.  We  will  not  show  the  distribution  of  reflected  and  transmitted  power 
since  our  analysis  will  proceed  after  replacing  QWR2  with  a  circular  retarder. 

Figure  25  shows  the  same  path  as  above  except  we  have  now  replaced  QWR3  with  a 
half  wave  circular  retarder,  CR2,  and  added  a  half  wave  linear  retarder.  We  replaced 
QWR3  with  CR2  for  the  benefit  of  beams  entering  module  1  from  the  signal  port  of  PBSl, 
and  we  will  explain  why  a  half  wave  retarder  has  been  added  next. 

CRl  rotates  the  linear  polarization  states  reflected  from  PBSl  45°,  and  the  plane 
mirror  before  the  ICU  flips  their  orientation  into  that  shown  in  Figure  25-V.  The  linear 
polarization  of  each  beam  is  oriented  at  45°  with  respect  to  the  s  and  p  directions  each 
beam  encounters  at  PBS2.  The  beams  intensity  of  each  beam  is  split  by  PBS2.  Half  the 
beam  is  transmitted  toward  the  plane  mirror  and  half  is  reflected  toward  the  roof  mirror. 
The  intensity  transmitted  toward  the  plane  mirror  is  plotted  as  a  function  of  object 
coordinate  in  Figure  26.  The  intensity  reflected  toward  the  roof  mirror  is  plotted  as  a 
function  of  object  coordinate  in  Figure  27.  The  Jones  calculus  for  this  path  is  provided  in 
Figure  28. 

As  stated  earlier,  the  path  from  PBS2  to  the  plane  mirror  and  back  to  PBS2  inside  the 
ICU  is  similar  to  the  path  analyzed  in  Figure  18,  and  its  performance  is  improved  when  a 
circular  retarder  is  employed. 

The  path  involving  the  roof  mirror  inside  the  ICU,  however,  benefits  from  the  use  of 
a  quarter  wave  retarder.  Figure  29  illustrates  the  ray  trace  involved  in  this  situation.  Since 
two  reflections  occur  with  the  roof  mirror  the  orientation  of  polarization  states  are  align 
with  the  p  states  of  PBS2  for  transmission  into  module  2.  The  Jones  calculus  for  this  path 
is  given  in  Figure  30. 

Figure  31  shows  the  path  from  the  PBS2  in  the  ICU  to  PBS3  in  module  2.  Note  that 
the  polarization  directions  emitted  from  PBS2  (Figure  31-IV)  are  not  aligned  with  the  p 
orientations  of  PBS3.  The  orientations  of  the  linear  polarization  of  these  beams  can  be 
aligned  with  the  p  directions  of  PBS3  by  inserting  a  half  wave  retarder  at  0°  anywhere 
between  PBS2  and  PBS3.  The  half  wave  retarder  flips  the  polarization  directions  about  the 
fast  axis  of  the  retarder  to  align  their  polarization  directions  with  the  p  directions  of  PBS3, 


Intensity  distribution  after 
transmission  through  PBS2 

(Plane  mirror  path  of  ICU;  design  utilizes  circular  retarders 

and  half  waveplates) 


*  Plotted  as  a  function  of  object  coordinate,  defined  at  PRl 
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Figure  26 


Jones  calculus  of  path 
from  s-SEED  to  ICU 

(utilizing  a  circular  retarder  and  a  halfwave  retarder) 
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Figure  30 


(Design  Altered) 
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III.  Polarization  IV.  Polarization  V.  Polarization  VI.  Polarization 

before  PBS2  after  PB.S.2  before  PBS3  after  PRS3 


The  remaining  paths,  through  PBS3  and  reflection  off  PBS3  are  similar  to  the 
situation  analyzed  in  Figure  18,  and  are  benefited  by  replacing  the  quarter  wave  linear 
retarders  with  half  wave  circular  retarders. 

The  changes  in  the  polarizing  element  design  suggested  by  this  analysis  are  shown  in 
Figure  23  and  can  be  cc’^trasted  with  the  existing  design  shown  in  Figure  13.  Figure  32 
shows  the  fraction  of  light  which  is  transmitted  by  the  entire  path  analyzed  of  the  existing 
Photonic  Switching  Network  as  a  function  of  object  coordinate  over  a  10°  FOV.  The  final 
result  was  obtained  by  summing  the  intensities  of  both  paths  in  the  ICU.  Figure  33  shows 
the  fraction  of  light  which  is  transmitted  by  the  entire  path  analyzed  of  the  modified 
Photonic  Switching  Network  as  a  function  of  object  coordinate  over  a  10°  FOV.  The 
improvement  of  efficiency  exists  only  for  off  axis  object  points  for  which  y5=0.  and  assumes 
ideal  circular  retarders  with  uniform  angular  response.  Research  is  needed  to  determine 
what  kind  of  circular  retarder  is  suitable  for  optical  interconnect  applications. 
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Figure  32 


3.  Mueller  matrix  imaging  polarimeter 

This  section  describes  the  Mueller  matrix  imaging  polarimeter.  After  a  summary  of 
Mueller  matrix  formalism,  the  Mueller  matrix  Imaging  Polarimeter,  operation,  data 
reduction  and  output  are  described. 

3.1  Mueller  Matrix  Formalism 

The  effect  of  an  optical  system  (element)  on  the  polarization  of  incident  light  is  most 
generally  expressed  by 

S  =MS  (1) 

where  S  and  S  are  the  4  x  1  Stokes  vectors  of  the  incident  light,  and  transmitted 
(reflected)  light,  and  M  is  the  4  x  4  real  Mueller  matrix  that  characterizes  the  linear 
interaction  of  the  light  with  an  optical  system.  The  elements  of  the  Stokes  vector  are 
defined  as 


(.  ^  p  .  ~  ^  )  \ 

(2:11  p.J:,  COSY) 
s.sin  Yly 


(2) 


where  E  ^  and  E  „  are  projections  of  the  electric  field  on  the  s  and  p  planes,  and  y  is  the 
difference  in  phase  accumulation  between  s  and  p  polarizations.  5o  is  the  total  intensity, 
i'  1  is  the  preference  for  p  polarized  light,  sa  is  the  preference  for  45°  polarization,  and  s  3 
is  the  preference  for  right  circularly  polarized  light.  The  inequality 


indicates  that  polarized  or  partially  polarized  light  may  be  described  by  a  Stokes  vector. 
The  polarization  of  a  beam  can  be  expressed  in  a  number  of  other  useful  ways  using  the 
Stokes  elements.  The  degree  of  polarization,  the  fraction  of  polarized  light  is  defined  as 


l)OP  = 


(-0 


I) OP  varies  from  zero  for  unpolarized  light  to  one  for  polarized  light.  Tlie  degree  of  linear 
polarization  is  the  fraction  of  linearly  polarized  light, 


DOLP  = 


:5) 
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DOL.P  is  zero  for  unpolarized  light  or  circularly  polarized  light  and  one  for  linearly 
polarized  light.  The  degree  of  circular  polarization  is  defined  by 

DOCP  =  —  (6) 

So, 


DOCP  is  zero  for  unpolarized  and  linearly  polarized  light  and  one  for  right  or  left 
circularly  polarized  light.  A  polarization  ellipse  may  be  drawn  from  the  eccentricity  and 
orientation.  The  eccentricity  is 


X  = 


I 


■arcsin 


V5t  +  S2  +  S3. 


(7) 


and  orientation  is 

1  ,  f 

a  =  -arctani  —  •. 

2  VS)' 


(8) 


For  a  non-depolarizing  optical  system,  the  Mueller  matrix  contains  seven  degrees  of 
freedom.  These  are: 

( t)  overall  intensity  transmission, 

(2,3)  linear  diattenuation  and  retardance  along  the  coordinate  axis. 

(4.5)  linear  diattenuation  and  retardance  at  45°  to  the  coordinate  axis, 

(6.7)  and  circular  diattenuation  and  retardance.  The  overall  phase  of  the  light  is  not 
preserved.  In  the  most.general  case,  the  Mueller  matrix  contains  16  degrees  of  freedom. 

The  Mueller  matrix  of  an  optical  system  (element)  is  a  function  of  wavelength  and 
angles  of  the  incident  light.  The  response  of  an  optical  system  to  a  spherical  wave  can  be 
described  by  a  Mueller  matrix  written  as  a  function  of  object  coordinate  /-/ ,  and  pupil 
coordinate  p  and  wavelength  X.  .12,13 


3.2  Mueller  matrix  Imaging  Polarimeter  Description 

The  Mueller  Matrix  Imaging  polarimeter  operates  by  sending  a  set  of  polarized 
spherical  waves  into  a  sample  and  measuring  the  polarization  state  across  the  transmitted 
or  reflected  spherical  wave.  It  takes  a  series  of  64  measurements  on  each  sample  and  uses 
Fourier  methods  to  reduce  the  data  to  provide  information  on  how  the  optical  system 
changes  the  polarization  of  light  on  a  ray  by  ray  basis.  From  the  data,  the  polarization  state 
of  a  ray  transmitted  from  a  sample  or  optical  system  can  be  determined  as  a  function  of  a 
ray's  incident  polarization  state  and  angle  of  incidence.  The  polarization  information  can 
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be  displayed  as  a  function  of  field  of  view  or  pupil  coordinate  in  the  form  of  a  Mueller 
Matrix.  Diattenuation  and  retardance  of  a  sample  or  optical  system  can  be  displayed  as  a 
function  of  field  of  view. 

The  Mueller  Matrix  Imaging  Polarimeter.  illustrated  in  Figure  34.  can  be  divided  into 
two  sections:  a  polarization  state  generator  and  a  polarization  state  analyzer.  The 
polarization  state  generator  consists  of  a  source,  a  He-Ne  633nm  laser  or  an  850nm  laser 
diode  and  polarizing  optics.  The  polarizing  optics  are  a  stationary  linear  polarizer  followed 
by  a  quarter  wave  linear  retarder  mounted  in  a  computer  controlled  rotary  stage. 

Following  the  polarizing  optics,  focusing  optics  prepare  a  spherical  wave  which  is  sent  to 
the  sample  under  test.  The  polarization  analyzer  consists  of  focusing  optics  followed  by  a 
quarter  wave  linear  retarder  followed  by  a  stationary  linear  polarizer  followed  by  a  ccd 
detector.  The  method  of  operation  is  as  follows.  The  first  linear  quarter  wave  retarder 
rotates  in  5.625°  increments  over  360°  and  the  second  quarter  wave  retarder  rotates  in 
increments  5  times  that  of  the  first,  28.125°.  The  two  rotate  together  sync.  The 
polarization  modulation  produced  by  the  discrete  rotations  of  the  retarders  of  the 
polarization  generator  and  analyzer  produces  an  intensity  modulation  measured  by  the  ccd 
detector  which  is  then  Fourier  analyzed.  In  this  case  all  16  elements  of  the  Mueller  matrix 
is  encoded  onto  the  detected  signal  at  each  pixel. 

The  theory  used  for  the  data  reduction  algorithm  of  the  Mueller  matrix  Imaging 
Polarimeter  is  given  next. 

Let  SC  / .  m  .  n)  =  (s'o.  .■?! .  s^.s  be  the  Stokes  vector  of  the  light  leaving  the 
polarization  generator.  \  =  in,  /  =  1.2.3.  1  be  the  Mueller  matrix  of  the 

sample  under  test,  and  11(C .  m.  n)  =  be  the  first  row  of  the  Mueller 

matrbc  of  the  polarization  analyzer.  The  direction  cosines  L.in.a  which  describe  the  ray 
directions,  will  be  suppressed  to  simplify  this  development.  The  detected  intensity  on  the 
ccd  detector  is  /  =  clliVI  S ,  where  cis  the  responsivity  of  the  ccd  detector.  Carrying  out  the 
matrix  multiplications  we  obtain 

•( 

/  =  c  I  (1) 

I.  /-  I 

where 

(1,.,  =  -'?..^',  (2) 
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When  the  rotational  step  of  the  polarization  analyzer  is  five  times  the  rotation  of  the 
polarization  generator,  [a  becomes 


[i-/< 


V 


cos^20 
sin  20'cos20' 
-  sin20 


cos  20 

cos*20cos^20 
cos^20sin20  cos20 
-cos^20sin  20' 


sin20cos20 
sin20cos20cos*20' 
sin  20cos20sin  20'cos20 
-sin20cos20sin20 


sin  20 

sin20cos^20' 
sin 20sin 20  cos20 
-  sin20sm  20 


(3) 

where  9  =  50  and  0  is  angle  between  the  fast  axis  of  the  quarter  wave  retarder  and  the 
transmission  axis  of  the  polarizer  of  the  polarization  generator  and  0  'is  the  angle  between 
the  fast  axis  of  the  polarization  analyzer  quarter  wave  retarder  and  the  transmission  axis  of 
the  polarization  generator  polarizer.  If  some  trigonometric  identities  are  applied  to  the 
result  in  Eqn.  3,  we  obtain  a  periodic  signal 

12 

/  =  ao+  Z  (a„cosn20  + {)„sinn20)  (4) 

n-  1 


Twenty  five  Fourier  amplitudes ao.an,b„,(n  =  1,2,... 12) are  determined  by 
performing  an  FFT  on  the  signal  I  to  obtain  the  Fourier  amplitudes  in  terms  of  the  scaled 
Mueller  matrix 

iVl'  =  ioiVl  (5) 

Lq-cIx  ' 

where  x  is  the  transmittance  of  the  polarizing  optics  and  focusing  optics  for  incident 
unpolarized  light. 

Note  that  four  Fourier  amplitudes  are  identically  zero,  and  that  five  identities  are 
satisfied  (0.7  =  -a3,b7  =  b3,a(,  =  a4,aii=-a<),b  ii=  -b^).  Deviations  from  these 
conditions  are  symptoms  of,  and  are  used  to  characterize  imperfections  in  the  polarimeter, 
such  as  azimuth-angle  errors,  non-ideal  components,  or  nonlinearity  of  detector  response. 
The  dependence  of  the  Fourier  coefficients  on  the  Mueller  matrix  elements  are 


48 


1  .  1  .  1  .  1- 
cto  =  2^ "  8 

a ,  =  0 

1-1- 
Ci2  =  -^fni2*Qin22 
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®2  13  ■^0^23 
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Ot  =  --m, 
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I  ■ 


1  . 
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1  .  I 


^8  "  1  ^^22' 


16  16 


^9  *  g^34 


10  “  ■*■  g'^22 


*^^11  "  ~3^^^34 


<^12“Yg^22  [5^33 
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^><5  =  ~:^/7i41  --m42 

O  O  1 


5^  =  0 

^7  “  ~  42 
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b<,  =  ~-mm24 
o 
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The  Fourier  coefficients  are  then  inverted  to  obtain  the  Mueller  matrix  elements 


m.^.i  =  2(,-a^-^aj 
77143  =  4(-a3  +  a6) 
^42  =  "^(^3-*-^7) 

I  •  1;, 

^41  ="2'^42-^i’s 


^33~^(*^8~*^^I2) 
^23  ~8(-bg4-5|2) 
77132  ~  S(6  g  4"  b  12) 

■  ,  I 

^12  =  '^^2- -77123 


/^24  =  "K-^’9-*'^1|) 

7n34  =  4(-a,-a,,) 


/71-  IT  “  •‘1’  T  rrt  ' 

10  ^  i 


49 


1  ■  I- 

/^i4  =  -2'^24"-46,--/n34 

/n22  =  8(a8-a,2) 


n\2\  ~  10  ~  2^22 

^31 


mna'i,  =  2ao--An,2--/n2i 

(7) 

To  obtain  the  absolute  Mueller  matrix  elements,  /rij.  ,  =  it  is  necessary  to 

measure  Lq.  This  is  normally  done  in  a  calibration  step  in  which  the  sample  is  removed  and 
the  polarimeter  is  operated  in  the  straight-through  configuration.  This  step  can  not  be 
done  in  an  imaging  polarimeter  for  most  samples.  Generally,  the  sample  introduces 
wavefront  aberrations  into  the  beam  causing  a  transverse  shift  in  ray  paths  when  it  is 
removed.  The  task  of  compensating  for  this  shift  in  software  would  be  formidable  and 
could  not  be  generalized  so  that  it  would  work  for  many  different  samples  with  different 
wavefront  aberrations.  We  instead  normalize  the  Mueller  matrix  by  dividing  all  the 
elements  by  m , ,.  Information  on  the  sample’s  transmittance  of  unpolarized  light  is 
therefore  lost. 

Accuracy  of  the  polarimeter  is  sensitive  function  of  the  alignment  errors  in  the 
quarter  wave  retarders  and  polarizers  and  the  retardance  errors  in  the  quarter  wave 
retarders.  The  alignment  errors  of  the  polarizing  optics  and  retardance  errors  in  the 
quarter  wave  retarders  can  be  measured  and  compensated  for  in  the  data  reduction.  The 
data  reduction  algorithms  of  the  Imaging  Polarimeter  include  algorithms  which 
compensate  for  these  errors.  A  single  channel  polarimeter  is  currently  being  constructed  to 
measure  the  alignment  errors  and  the  retardance  errors  of  the  polarizing  optics  of  the 
Imaging  Polarimeter.  The  full  analysis  and  the  single  channel  polarimeter  is  described  in 
section  5. 


4.  Stokes  Imaging  Poiarimeter 

Imaging  Polarimetry  supplements  phase  measuring  interferometry  by  measuring  the 
polarization  content  over  a  wavefront.  A  Stokes  Imaging  Poiarimeter  measures  Stokes 
vectors  over  the  wavefront.  The  UAH  Imaging  Poiarimeter  is  operated,  in  one 
configuration,  as  a  Stokes  Imaging  Poiarimeter  In  this  configuration  we  introduce  a 
plane  wave  or  spherical  wave  of  known  polarization  state  to  the  optical  element  or  optical 
system  under  test,  and  measure  the  Stokes  vectors  of  the  transmitted  light  on  a  ray  by  ray 
basis.  We  characterize  the  response  of  the  optical  system  to  a  waveform  of  known 
polarization. 

The  Stokes  Imaging  Poiarimeter  presently  operates  in  the  visible  and  near  infra-red. 
The  instrument  has  been  set  up  and  is  operated  in  the  Research  Institute  building  at  the 
University  of  Alabama  in  Huntsville.  The  instrument  is  fully  automated  with  computer 
controlled  rotary  stages  and  involves  polarizers,  retarders,  lenses  and  a  ccd  camera.  Once  a 
data  collection  run  is  started  by  the  controlling  computer,  the  automated  stages  are  rotated 
in  increments  according  to  a  rotating  retarder  scheme  described  below  and  several  images 
are  digitized  using  a  image  frame  grabbing  board.  At  the  end  of  the  data  collection  routine 
the  modulated  intensity  images  are  shipped  over  the  PC  bus  and  Fourier  analyzed  to  obtain 
the  Stokes  elements  in  image  form.  The  Stokes  images  are  then  manipulated  to  obtain 
relevant  polarization  images. 

Stokes  vectors  across  a  wavefront  or  ’Stokes  Images’  are  generally  measured  at  the 
exit  pupil  by  imaging  the  exit  pupil  the  ccd  detector.  When  the  exit  pupil  is  imaged  onto 
the  ccd,  each  pbcel  of  the  ccd  gathers  light  from  a  small  bundle  of  rays  that  have  passed 
through  a  sample  and  the  Imaging  Poiarimeter  measures  the  polarization  aberration 
functionl7»lS.  An  image  plane  of  an  optical  system  can  also  be  imaged  onto  the  ccd, 
wherein  the  distribution  of  polarization  states  in  the  point  spread  function  is  measured  and 
the  polarization  point  spread  function  of  the  system  is  determined  in  Stokes  vector  form. 

The  objective  of  the  imaging  poiarimeter  is  to  determine  a  Stokes  vector  associated 
with  each  pixel  of  an  input  scene  or  wavefront.  The  configuration  of  the  poiarimeter  is 
given  in  Figure  35.  It  consists  of  a  rotating  retarder  poiarimeter  followed  by  a  ccd  camera 
connected  to  an  image  processor.  The  poiarimeter  contains  a  quarter  wave  retarder 
(QWR)  and  a  linear  polarizer.  The  poiarimeter  operates  by  stepping  the  quarter  wave 
retarder  through  a  360°  cycle  in  11.25°  increments  for  a  total  of  32  sequentially  measured 
frames.  The  resolution  of  each  measured  frame  is  reduced  by  averaging  5  by  4  pbcel 
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Figure  35  Schematic  diagram  of  UAH  Stokes  Imaging  Polarimeter. 


squares  into  bixels.  The  averaging  is  done  to  increase  signal  to  noise  and  to  minimize  the 
effect  of  beam  wander.  The  resolution  of  each  measured  frame  is  reduced  after  each 
measurement  to  minimize  memory  requirements.  Beam  wander  is  discussed  in  the  error 
and  calibration  section.  Producing  a  set  of  bbcels  over  the  measurement  region.  A  Fast 
Fourier  transform  is  performed  on  the  data  bixel  by  bixel  to  generate  the  Fourier 
coefficients  associated  with  each  bixel  of  the  input  image.  The  Stokes  parameters  can  then 
be  obtained  from  the  Fourier  coefficients  on  a  bixel  by  bbcel  basis.  The  Fourier  analysis 
related  to  the  Stokes  parameters  is  described  below. 

The  Mueller  matrix  of  an  ideal  quarter  wave  retarder  and  ideal  polarizer 
combination  is 


M(0)  =  - 


1  cos‘2G  cos20sin20  -sin20 
1  cos^20  cos20sin20  -sin20 


where  0  is  the  angle  between  the  fast  axis  of  the  quarter  wave  retarder  and  the  transmission 
axis  of  the  polarizer.  Since  the  polarizer  remains  stationary  before  the  ccd,  only  one  state 
of  linear  polarization  is  incident  on  the  ccd.  The  polarization  sensitivity  of  the  ccd  detector 
does  not  have  to  be  considered. 

The  detector  measures  intensity,  /.  The  intensity  of  the  light  imaged  onto  the 
detector  is  modulated  by  the  rotating  retarder  according  to 


/(O)  =  So  -pC^'o  +  s,  cos'" 20  + S2 cos 20 sin 20 -S3 sin  20).  (2) 

or  rewriting  using  trigonometric  substitutions 

/(0)  =  So  =  ^(^So  +  ^s,  -S3sin20  +  ^s,cos40  +  ^S2sm40j,  (3) 

where  0  is  the  rotation  angle  of  the  quarter  wave  retarder  relative  to  the  transmission  axis 
of  the  polarizer.  Eqn.  3  is  a  Fourier  series  with  Fourier  amplitudes 


^^0 “ 2V  2^' 

t>2  —  ~S3 , 


“4= 0^1 


S3 


(4) 


The  Fourier  coefficients  are  inverted  to  obtain  the  Stokes  vector  elements 

S Q  =  2.CLq~  CL 

=  2a. ^ 

.92  =  26, 

S3  =  -62.  (5) 

The  intensity  across  the  scene  or  wavefront  is  divided  out  by  dividing  Sj, 93,93  by  So 
bixel  by  bixel.  The  degree  of  linear  polarization  (DOP),  degree  of  linear  polarization 
(DOLP),  degree  of  circular  polarization  (DOCP)  can  also  be  displayed  as  a  function  of 
field  of  view.  The  eccentricity  and  orientation  of  a  polarization  ellipse  can  also  be 
displayed.  They  are  a  series  of  polarization  ellipses  drawn  over  the  pupil. 

In  another  configuration,  the  Stokes  imaging  polarimeter  uses  a  rotating  polarizer  in 
front  of  the  ccd  camera.  Three  images  are  taken  with  the  Stokes  polarimeter  in  this 
configuration  and  the  data  is  reduced  to  give  the  first  three  stokes  parameters:  9o ,  9 1 , 93. 
The  usefulness  of  the  Stokes  Imaging  Polarimeter  in  this  configuration  is  to  take  the  data 
very  quickly,  2.5  seconds.  It  is  used  to  take  polarimetric  images  of  non-static  scenes.  The 
Stokes  images  can  be  manipulated  to  obtain  all  of  the  above  mentioned  polarization 
quantities,  DOP,  DOLP,  etc. 
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5.  Single  Channel  Mueller  Matrix  Polarimeter 

A  single  channel  Mueller  matrix  polarimeter  is  currently  being  constructed  in  the 
Research  Institute  building  at  the  University  of  Alabama  in  Huntsville.  A  schematic  of  the 
polarimeter  is  shown  in  Figure  36.  The  instrument  will  support  our  research  in  Imaging 
Polarimetry  by  1.  calibration  of  polarizing  optics  used  in  the  imaging  polarimeter,  2. 
measuring  alignment  errors  of  the  quarter  wave  retarders  and  polarizers  in  the  Imaging 
Polarimeter,  3.  Accurate  polarization  measurements  of  interesting,  4.  Depolarization 
studies  of  thin  films,  psuedo-depolarizers,  and  other  scientifically  valuable  samples. 

The  instrument  is  fully  automated  with  computer  controlled  rotary  stages  and 
involves  polarizers,  retarders,  lenses  and  a  silicon  photodiode  broadband  photo-detector. 
The  photo-diode  is  connected  to  a  pre-amp  connected  to  a  4  KHz  12-  bit  A-D  converter 
board.  Once  a  data  collection  run  is  started  by  the  controlling  computer,  the  automated 
stages  are  rotated  in  increments  according  to  the  same  rotating  retarder  scheme  as  is  used 
in  the  Irnaging  Polarimeter  and  several  data  points  are  collected.  At  the  end  of  the  data 
collection  routine  the  modulated  intensity  data  points  are  shipped  over  the  PC  bus  and 
Fourier  analyzed.  In  one  configuration,  the  Single  Channel  Polarimeter  is  used  to  calibrate 
the  polarizing  optics  of  the  Imaging  Polarimeter.  The  calibration  run  gives  orientation 
errors  of  the  retarders  and  polarizers  of  the  Imaging  Polarimeter  and  retardance  errors  of 
the  retarders.  In  another  configuration,  it  has  the  capacity  to  accurately  measure  the 
Mueller  matrix  of  a  sample. 

5.1  Polarizing  element  orientation  and  retardance  error  compensation 

The  procedure  to  find  orientation  errors  and  retardance  errors  of  the  Imaging 
Polarimeter  polarizing  optics  will  be  described  first.  The  errors  discussed  here  are  the 
retardance  errors  of  the  first  and  second  retarders,  e  i  and  62.  orientation  of  the  first  and 
second  quarter  wave  retarders,  63  and  €4,  and  the  orientation  error  of  the  last  polarizer  65. 
The  orientation  errors,  63 ,  e4 ,  G5  are  referenced  to  the  transmission  axis  of  the  first 
polarizer.  The  five  errors  are  shown  graphically  in  Figure  37. 

Calibration  of  the  Imaging  Polarimeter  is  extremely  important  since  orientation 
alignment  errors  of  the  polarizing  optics  introduce  error  into  the  measured  Mueller 
matrices.  Misalignment  of  the  first  retarder  relative  to  the  transmission  axis  of  the  first 
polarizer  by  1°  can  cause  an  error  of  up  to  7%;  misalignment  of  the  second  retarder  by  the 
same  amount  causes  up  to  34%  error  and  retardance  error  in  either  of  the  two  retarders 
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TA  =  transmission  axis  of  linear  polarizer 


produces  1%  error.  []  The  orientation  and  retardance  errors  e,,  i  =  1,2, 3, 4, 5  can  be 
included  into  the  equations  relating  the  Fourier  amplitudes  and  the  Mueller  matrix 
elements.  The  following  notation  will  be  used: 


a,  -2e4-2e3-26g 
a2  =  2e4  +  2e3-2es 
*^3~2e4-4e3-26s 
a4  =  2e4  +  4e3-2e5 
^5  ~  265  —  264 
^6~265-464  +  263 

3i  =  1  -  COS 6, 

(Bg  =  1  -  COS62 

6,.2  =  90‘’  +  e,.2 


0t7~26c~4£4-'26- 


otfl  =  -a, 


a9  =  4e4-4e3-265 
aio  =  4e4  +  463-2e5 
^11  ~  464- 265 


^3  =  1  +  cos6, 
(34  =  1  +  COS62 


The  following  equations  relate  the  Fourier  amplitudes  to  the  Mueller  matrix 
elements  and  errors: 

ao=^;m,,^;p3m,2  +  J(34Cos26sm2,-4-i|33p,cos265m22  +  i34sin2esm3,  +  i33p,sin26sm32 

a,«isin6,sin2e3m,4  +  jsin6,134Sin,2€3.cos2es/n24  +  4-34Sin6,sin2e3Sin265/n34 
«,-rP.«S‘'*,«.:-;P.s.n4«,m„-ip.0.ccs4.,cos2c,m„*lp,p.sin4c,cos2.,m„.lp,p,co.4c,sin2c,m„.lp,p.sm4c,s.n2e,m„ 

^3  8PiSin62sina3m42-g(3i  sin 62Cosa3m43 
a4  =  --sin6,sin62Cosaim4.^5  =  isin62Sinasm4,4.i(33sin62Sina5m42 
a6  =  -mn6,sin62COsa2/n44a7  =  -i(3,sin62sina4m42  +  5f3,sin62Cosa4/n43 
=  IiPiP2Cosa9(m22  +  /n33)4-^(3jP3Sina9(m32-/n23) 

*^9  =  5|32sin6isina6/n24  +  g(32sin6,  cos  77134 

aio  =  4|B2Cosa,,m2, +||32(33Cosa,,m22  +  jl32Sina,,m3, +Ip233sina,,n732 

^11  =  ~  iPaSin  6  j  sin  a7  77124 -5(32  sin  6  jcostty  77134 

<^12  =  T6PlP2COSa,o(77l22  -  77X33)  +  ^f3jp2Sina,o(77l23  +  /7l32) 
h,=5sin6,cos2£377z,,4  +  i(3,sin6,cos2e3cos2€577i24-^i|34sin6,cos2e3s]n2es77i34 
t.:.-lP.s.n4c,m.,.lp.cos4c,m.,.lp,p.cos4e,cos2c,m„-lp,p.si„4c,cos2.,m„.lp,p.cos4c,s,„2c,^„-ip.p.„„.,c,s.n2c,m„ 

^3  “  “sPl  S2COSaj77l42  +  Sin  62Sina3  77l43 

64  =  -sin6,sin62sina,77i44h5  =  -jsin62Cosa577i4,-i|S3sin62Cosa577i42 
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be  “  “^sin6,  sin62sina2m44b7  “gPi  sin62Cosa4/n42-5Pi  sin62sina4^43 

i>8  =  -r5l3ip2Sina9(m22  +  /n33)-^p,(32C0sa9(m23-m32) 

t»9  =  -  5|32sin  6,  cosae/7i24  +  gP2sin  6,  sinae/n34 

i’  io  =  -5P2sina,ifn2i  -5P2P3Sina,,/n22  +  ^p2C0sanna3,  +  ^(32|33Cosai,m32 

t>ii  =  5P2sin6|  cosa7m24-g32Sin6, sina7m34 

f>i2  =  -i7PiP2sina,o(na22-^33)  +  T^I3i(32COsa,o(m32  +  m23) 


which  can  then  be  inverted  to  obtain  the  Mueller  matrix  elements  in  terms  of  the 
Fourier  amplitudes  and  errors 


_  8 

^44  3sin 6 ,  sin  67 


cosai  cosa^  sin 


i;)”43=^ 


-yl3COsa3''S3Sina3*/<7COsa^-S7Sin 
8 1  sin  67 


fUAy  =  -8 


^3Sina3-S3C0SO3->/l7sino^''S7C0sa^  ’♦fl; 


(S I  sin  67 


2  cos  05  sin  67 


/i^cos Oj^-S^sin  Oj'  /I  II  cosa^-fin  sin  07 


|}2Sin  6| 


^2  sin  6 1 


-P4Cos2es'n24  .4B|  P4sin2€5m34 

/ni4  =  ^  cos 263 sin 8 1  ^~2  ^22  “16 


/i  g  cos  a, -/( 12  cos  a  IQ- 0g  sin  a,  *5 12  sin  a, 9  - /Igsin  a,* /i  |2Sin  o  ig-BgCosa,- B  |2C0sa  |q 

^33=  16 - ^ - — - ^23=  16 - - 


/n37  =  16 


/igSin<i2-/l|2Sina  ,q-  BgCosa,-B  |2C03a  |q 


in,7  = 


16/12 cos 4C3-I6S2 sin  4€3-p  I P4 cos 2€5m22"P I P4 sin 265(7132 


16 /1 2  sin  463-  16S2<:os  463-P  1P4  cos  2  £5(71 33- P1P4  sin  16/1  ujcosa  n  -  16B  i^sin  a  n  -p2p3/n22 

/71i3  ^  2?I  ^21  2P^  ■ 


ni3l  = 


-(P2P3'"32-  16B|oCosa||  -  16/1  igsin  Oil) 


m,,  =4/lo-^P3^i2“iP4COs26sm2i  -JI33P4COS 2650122-5134 sin 265^31  -  j|33p4Sin 2esOZ32 

(7) 
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where  the  Mueller  matrix  elements  should  be  calculated  in  the  order  shown.  These 
equations  must  be  solved  sequentially  because  some  of  the  matrix  elements  have  been  left 
in  terms  of  other  matrix  elements  to  reduce  the  complexity  of  the  equations. 

To  solve  for  the  five  alignment  and  retardance  errors,  the  single  channel  Mueller 
matrix  polarimeter  is  operated  with  no  sample  and  the  Fourier  amplitudes  are  obtained 
from  the  measured,  modulated  intensity.  Elements  of  the  identity  matrix  are  inserted  for 
the  Mueller  matrix  elements  in  Eqns  7,  yielding 

ao  =  i  +  5p3p4Cos2€5 

az  =  5|3,(34COs4€3Cos2es  +  5(3,(34sin4€3Sin2e^4  =  -^sin6,  sin  63  cos  a, 

06  =  ^ sin  6 1  sin  62 cos Oj 

Cts  ~  g  P 1  P2  *^9 
10  ~  5P2P3  cos  a ,  1 
a, 2  =  0 


.  62  =  -gPiP4Sin4e3Cos2e5  +  g(3,|34Cos4e3sin 26^4  =  ^sin6,sin62Sina, 
b6  =  ~5sin6isin62sina2 

6g  =  -5P,P2Sina9 

b  10  “  "gPsPsSinaii 

6 12  =  0 

(8) 


Note  that  all  odd  Fourier  coefficients  are  set  to  zero,  since  their  presence  is  possible  only 
when  alignment  errors  are  present.  These  equations  are  then  inverted  to  obtain  the 
alignment  and  retardance  errors. 


Gs  =  -tan' 


-  .i/'aiocosa,-a8cosa,,’\ 

61= COS  - 

*  \  a  locosq^^agcosaji  J 
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(9) 


62  =  cos 


->r 

V 


•*°  10 
COS  a  II 


Once  these  equations  are  obtained,  they  are  inserted  back  into  Eqns  7  to  compute  the  error 
corrected  Mueller  matrix  elements. 
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6.  Measurements 

This  section  presents  several  measurements  made  with  the  imaging  polarimeter.  The 
purpose  of  this  section  is  to  convey  the  capabilities  of  the  imaging  polarimeter.  The 
measured  samples  and  results  will  only  be  briefly  discussed.  The  measurements  made  on 
the  polarizing  beam  splitter  cubes  are  presented- in  detail  in  Appendix  B,  and  will  not  be 
shown  in  this  section. 

Before  each  measurement  cycle  the  single  channel  polarimeter  is  operated  with  no 
sample  in  the  compartment  to  determine  the  initial  alignment  errors  of  the  quarterwave 
retarders  and  the  polarization  analyzer  linear  polarizer.  This  is  done  by  inserting  the  single 
charmel  polarimeter  detector  in  an  image  plane  of  the  imaging  polarimeter  after  the 
polarization  analyzer.  The  transmission  axis  of  the  polarization  generator  linear  polarizer 
defines  the  axis  in  which  orientation  errors  of  the  retarders  and  the  last  polarizer  are 
measured.  Also  in  this  step  we  check  for  retardance  errors  of  the  quarterwave  retarders. 
The  orientation  errors  and  retardance  errors  can  then  be  compensated  for  in  the  reduction 
of  data  from  the  imaging  polarimeter. 

A  calibration  Mueller  matrix  taken  before  a  data  run  is  shown  in  Figure  38,  The 
imaging  polarimeter  was  operating  over  a  =  -6®  field  of  view.  Since  no  sample  is  present  in 
the  sample  compartment  during  the  calibration  run,  the  Identity  matrix  is  expected,  i.e,,  all 
diagonal  elements  are  one  and  all  off  diagonal  elements  are  zero.  The  measured 
calibration  matrix  falls  within  the  expected  values  within  ±2,5% . 

The  Mueller  matrbc  of  a  total  internal  reflection  linear  polarizer  frf'.h  Optical 
Research  Associates,  measured  by  the  imaging  polarimeter,  measured  at  633  nm  is  shown 
in  Figure  39.  The  diattenuation  of  the  polarizer  can  be  extracted  from  the  9  upper  left 
Mueller  matrix  elements.  The  diattenuation  will  not  be  presented  here.  The  Mueller 
matrbc  was  measured  over  a  ±5°  field  of  view.  The  diattenuation  remains  high  up  to  3°, 
where  it  falls  off  dramatically,  which  is  marked  by  the  sharp  drop  off  in  one  direction  on  the 
four  elements  in  the  upper  left  corner  of  the  matrbc. 

Figure  40  shows  the  Mueller  matrbc  of  a  Meadowlark  Inc.,  liquid  crystal  circular 
polarizer,  measured  at  850nm.  The  accuracy  of  this  matrix  is  drastically  reduced  because 
the  measurement  technique  of  the  imaging  polarimeter  does  not  lend  itself  to  measuring 
circular  diattenuators.  The  reason  is  that  the  polarization  generator  sends  a  series  of 
polarization  states  that  are  orthogonal  eigenpolarization  states  of  the  circular  polarizer, 
limiting  the  amount  of  light  passed  by  the  sample.  The  loss  of  signal  for  these 
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measurements  introduces  substantial  amounts  of  noise  into  the  measurement.  The 
behavior  of  the  circular  polarizer,  however  is  present  in  the  lower  right  and  lower  left 
elements  of  the  matrix.  Further  research  is  needed  to  modify  the  present  measurement 
technique  to  allow  the  imaging  polarimeter  to  characterize,  circular,  diattenuators. 

A  measurement  of  the  uniformity  of  a  zero  order  linear  quarterwave  retarder  is 
shown  in  Figure  41.  The  Stokes  imaging  polarimeter  was  used  to  measure  this  sample  at 
633nm.  A  plane  wave  of  uniform  linear  polarization  was  incident  on  the  sample  such  that 
the  linear  polarization  direction  was  aligned  with  the  fast  axis  of  the  quarterwave  retarder. 
The  polarization  states  of  the  transmitted  co-linear  rays  are  plotted  as  a  function  of 
wavefront  coordinate.  The  transmitted  wavefront  polarization  is  nearly  uniform.  The 
non-uniformity  is  primarily  due  to  a  spatial  variation  of  the  crystal  axis  orientation  across 
the  sample.  Thickness  variations  also  contribute  to  non-uniformity  although  the  two  effects 
can  not  be  decoupled. 

7.  Calibration  Issues 

In  this- section,  the  issues  of  camera  calibration,  bootstrap  calibration  of  the 
polarization  generator  and  analyzer  elements,  movement  in  the  scene,  source/detection 
fluctuations,  beam  wander,  and  instrumental  polarization  are  discussed. 

7.1  Camera  calibration 

Camera  calibration  can  be  relaxed  somewhat  due  to  the  nature  of  our  measurement 
technique.  Our  calculated  images,  or  pupil  maps  in  the  case  of  the  PBSC  measurements, 
consist  of  added  and  subtracted  images  divided  by  another  image.  Normally  each  pixel 
grey  value  would  need  to  multiplied  by  a  sensitivity  factor  which  is  slightly  different  for 
each  pbcel  (sometimes  10%  difference),  but  since  our  images  ratioed,  the  sensitivity  factor 
divides  out.  However,  the  dark  count  (pixel  value  at  zero  light  level)  must  be  subtracted 
out  of  each  measurement.  All  data  reduction  algorithms  do  this. 

7.2  Polarizing  element  calibration 

Without  sufficient  prior  knowledge  of  the  properties  of  the  polarization  elements,  a 
bootstrap  calibration  needs  to  performed,  i.e.  a  calibration  using  retarders  and  polarizers 
with  unknown  properties.  The  calibration  should  make  the  fewest  number  of  assumptions 
about  the  polarizing  elements  and  utilize  methods  that  are  relatively  insensitive  to  the 
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Polarization  map 


Polarization  map  of  wavefront  emitted  from  a  quartz  quarterwave  retarder  as  a 
function  of  wavefront  coordinate.  Retarder  illuminated  with  a  plane  wave  of  linear 
polarization  parallel  to  the  fast  axis. 
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weakest  assumptions.  Once  an  accurate  Mueller  matrix  polarimeter  is  operating,  it  can  be 
used  to  calibrate  polarization  elements  for  other  polarimeters,  greatly  simplifying  the 
exacting  task  of  construction  and  accurate  imaging  polarimeter. 

Preliminary  calibration  has  been  done  on  the  UAH  imaging  Stokes  polarimeter. 
Polarizer  calibration  assumed  that  the  polarizers  only  display  diattenuation.  The 
transmittance  of  the  transmission  axis  and  the  absorption  axis  of  polarizers  were  measured. 
Two  similar  polarizers  were  calibrated  in  this  way  were  then  used  to  calibrate  the  retarders. 
The  first  nine  elements  m  i ,  through  of  the  retarder  Mueller  matrix  were  measured  by 
rotating  the  polarizers  at  a  3:1  ratio  on  either  side  of  the  quarter  wave  retarder.  The 
resulting  intensity  modulation  has  four  harmonics  and  a  dc  term  in  its  Fourier  series  which 
will  solve  for  the  nine  Mueller  matrix  elements.  The  linear  retardance  and  linear 
diattenuation  can  be  obtained  from  these  nine  Mueller  matrix  components  assuming  there 
is  no  circular  retardance.  Circular  retardance  occurs  in  multiple  element  linear  retarders 
when  the  fast  axes  of  the  linearly  birefringent  elements  are  not  exactly  0°  or  90°  apart.  To 
test  for  circular  retardance,  the  QWLR  was  placed  between  two  crossed  linear  polarizers 
with  its  fast  parallel  to  the  transmission  axis  of  the  first  polarizer.  If  circular  retardance  is 
present,  the  electric  field  vector  is  rotated  slightly,  so  that  the  final  polarizer  must  be 
rotated  to  locate  the  null.  No  circular  retardance  was  observed  in  these  retarders,  but  we 
have  detected  it  in  other  linear  retarders.  If  circular  retardance  is  observed,  the  remaining 
seven  components  of  the  QWLR  Mueller  matrix  should  be  measured.  A  method  to 
measure  these  components  is  described  in  reference  22. 

7.3  Speed  of  data  acquisition 

The  imaging  polarimeters  discussed  here  take  measurements  sequentially  in  time  and 
are  sensitive  to  drift.  Any  change  in  source  intensity  or  any  drift  in  camera  sensitivity 
during  the  measurement  sequence  is  a  source  of  1/f  noise.  The  measurements  must  be 
made  as  quickly  as  mechanically  possible.  The  rotation  stage  which  rotates  the  polarization 
generating  and  analyzing  elements  is  the  limiting  factor.  Ideally  the  stage  should  rotate  the 
polarizing  element  into  position  during  one  frame  period  or  l/30th  second.  This  would  be 
a  dead  frame  period.  In  the  case  of  the  Stokes  polarimeter,  if  frame  averaging  is  not  done, 
the  measurement  time  required  would  be  four  frames  for  data  acquisition  and  three  dead 
frames,  or  7/30th  of  a  second.  High  speed  data  acquisition  may  enable  polarimetric 
measurements  of  quasi-still  life  such  as  planes,  boats  or  vegetation. 
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7.4  Beam  wander 

Another  important  calibration  issue  is  beam  wander.  To  understand  the  problem, 
consider  how  digital  image  processors  calculate  derivatives.  A  partial  derivative  of  the 
image  intensity  in  the  x  direction  is  calculated  by  shifting  the  image  in  the  x  direction  and 
subtracting  the  image  from  itself.  The  imaging  polarimeter  performs  data  reduction  by 
adding  and  subtracting  images.  A  mo*,  'ng  object  in  the  scene  introduces  spurious 
contributions  into  the  Stokes  images  proportional  to  the  resulting  intensity  derivatives. 
These  intensity  gradients  are  often  greatest  around  the  objects  edges,  but  edges  are  where 
many  objects  have  their  greatest ;  olarization  signatures.  Fortunately,  in  optical  system 
polarimetric  metrology,  our  subject  is  usually  a  fixed  optical  system  and  the  intensity 
gradients  are  generally  small.  Thus,  far  more  accurate  polarimetry  is  possible  in  this  setting 
than  in  remote  sensing  or  tactical  polarimetry. 

The. faces  of  'QWLR  need  to  be  highly  parallel  because  wedge  in  a  rotating  element 
generates  beam  wander.  If  the  retarder  is  slightly  prismatic,  the  image  will  move  in  a  circle 
on  the  ccd  resulting  in  a  spatial  misregistration  and  imago  derivatives.  Large  amounts  of 
beam  wander  can  be  partially  compensated  by  carefully  shifting  the  measured  images  in 
software. 

7.5  Instrumental  polarization 

The  instrumental  polarization  due  to  wavefront  shaping  optics  within  the  polarimeter 
may  require  compensation  during  data  reduction.  Their  instrumental  polarization  will 
have  to  be  determined  each  time  a  new  optical  system  is  configured  within  the  sample 
compartment,  so  standardized  methods  for  removing  this  systematic  error  should  be 
developed. 
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8.  Conclusions 

This  research  has  produced  analytical  methods  to  aid  in  the  polarization  design  of 
polarization  based  optical  computers  and  has  developed  needed  instrumentation  to 
measure  and  understand  the  propagation  of  polarized  light  in  polarization  based  optical 
interconnects. 

We  provide  a  summary  of  the  polarization  analysis  done  on  the  AT&T  Iteration  3 
Photonic  Switching  Network.  The  analytical  approach  can  be  used  to  analyze  other 
polarization  based  optical  computing  architectures.  In  the  course  of  the  analysis  and 
research,  we  discovered  several  methods  which  may  improve  the  efficiency  of  the  optical 
interconnects  associated  with  AT&T  optical  computers.  The  changes  to  the  existing 
designs  involve  replacing  several  quarter  wave  linear  retarders  with  circular  retarders,  and 
inserting  half  wave  retarder  plates  in  two  places  to  manipulate  the  instrumental 
polarization  so  that  the  instrumental  polarization  works  in  our  favor. 

An  Imaging  Polarimeter  has  been  built  and  calibrated  to  perform  polarization 
metrology  on  polarizing  optical  systems.  It  presently  operates  at  633nm  and  849nm.  The 
Imaging  Polarimeter  can  be  operated  in  two  configurations.  In  the  first  configuration,  it 
measures  polarization  content  over  a  wave  front.  The  polarization  of  an  optical  system  can 
be  measured  at  the  exit  pupil  to  obtain  the  polarization  aberration,  or  it  can  measure  the 
polarization  content  of  anjmage  to  obtain  the  polarization  point  spread  function.  In  the 
second  capacity  it  sends  in  polarization  states  and  measures  output  polarization  states  to 
measure  the  polarization  properties  of  the  optical  system  for  a  set  of  ray  paths.  From  a  set 
of  measurements,  the  Mueller  matrix  of  an  optical  system  can  be  obtained  as  a  function  of 
field  of  view  and  pupil  coordinate. 

The  Imaging  polarimeter  demonstrated  its  usefulness  in  a  set  of  measurements  made 
on  polarizing  beam  splitter  cubes.  The  paper  which  developed  out  of  this  work  "Angular 
dependence  of  polarizing  beam  splitter  cubes,"  has  been  submitted  to  Applied  Optics  and  is 
attached  to  the  end  of  the  report.  Several  other  measurements  have  been  made  with  the 
Mueller  matrk  Imaging  Polarimeter,  such  as  quarter  wave  retarders,  circular  polarizers, 
and  total  internal  reflection  polarizers. 

Further  research  is  needed  to  increase  the  accuracy  of  the  Imaging  Polarimeter. 
Presently,  accuracy  of  the  Imaging  Polarimeter  is  limited  by  the  8-bit  A-D  conversion  of  the 
video  frame  grabber  board.  An  upgrade  to  a  12-bit  A-D  conversion  could  increase  the 
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accuracy  by  as  much  as  a  factor  of  10.  The  added  accuracy  is  needed  to  study  smaller 
polarization  effects,  such  as  weak  instrumental  polarization,  or  weak  field  of  view 
dependence  of  a  polarizing  device,  or  even  depolarization. 


7.  References 


1.  A.  Huange,  "Architectural  considerations  involved  in  the  design  of  an  optical  digital 
computer,"  Proc.  IEEE  72, 780-786  (1984). 

2.  D.A.B.  Miller,  D.S.  Chemla,  T.C.  Damen,  A.C.  Gossard,  W.  Wiegmann,  T.H.  Wood 
and  C.A.  Burrus,  "Novel  hybrid  optically  bistable  switch:  The  quantum  well 
selfrelector-optic  effect  device,"  Appl.  Phys.  Lett.  45, 13  (1984). 

3.  D.A.B.  Miller,  D.S.  Chemla,  T.C.  Damen,  T.H.  Wood,  C.A.  Burrus,  A.C.  Gossard  and 
W.  Wiegmann,  "The  quantum  well  self-electro-optic-effect  device , optoelectronic  bistability 
and  oscillation,  and  self  linearized  modulation,"  IEEE  J.  Quantum  Electron,  QE-21,1462 
(1985). 

4.  D.A.B.  Miller,  J.E.  Henry,  A.C.  Gossard,  and  J.H.  English,  "Integrated  quantum  well 
self-electro-optic  effect  device:  2x2  array  of  optically  bistable  switches,"  Appl.  Phys.  Lett. 

49, 821  (1986). 

5.  G.  Livescu,  D.A.B.  Miller,  J.E.  Henry,  A.C.  Gossard,  and  J.H.  English,  "Spatial  light 
modulator  and  optical  dynamic  memory  using  integrated  symmetric  self-electro-optic  effect 
devices,"  in  Postdeadline  Papers,  Conference  on  Lasers  and  Electro-optics  (Optical  Society  of 
America,  Washington  DC  1987),  paper  ThUll;  Optics  Lett.,  April  1988. 

6.  D.A.B.  Miller,  D.S.  Chemla,  T.C.  Damen,  A.C.  Gossard,  W.  Wiegmann,  T.H.  Wood 
and  C.A.  Burrus,  "Band-edge  electroabsorption  in  quantum  well  structures:  The  quantum 
confined  stark  effect,"  Phys.  Rev.  Lett.  53, 2173  (1984). 

7.  D.A.B.  Miller,  D.S.  Chemla,  T.C.  Damen,  A.C.  Gossard,  W.  Wiegmann,  T.H.  Wood 
and  C.A.  Burrus,  "Electric  field  dependence  of  optical  absorption  near  the  band  gap  of 
quantum  well  structures,  "Phys.  Rev.  B  32, 1043  (1985). 


72 


8.  A.L.  Lentine,  H.S.  Hinton,  D.A.B.  Miller,  J.E.  Henry,  J.E.  Cunningham,  and  L.M.F. 
Chirovsky,  "Symmetric  self-electro-optic  effect  devicerOptical  set-reset  latch, "Appl.  Phys. 
Lett.,  52, 1419  (1988). 

9.  F.B.  McCormick,  A.L.  Lentine,  L.M.F.  Chirovsky,  and  LA.  D’Asaro,  "An  all-optical 
register  using  symmetric  self  electro-optic  effect  devices,"  in  Conference  on  Photonic 
Switching  (Optical  Society  of  America,  Salt  Lake  City,  Utah,  1989),  paper  ThC5-l. 

10.  A.L  Lentine,  H.S.  Hinton,  D.A.B.  Miller,  J.E.  Henry,  J.E.  Cunningham,  and  L.M.F. 
Chirovskiy,  "Symmetric  self-elctro-optic  effect  device:  Optical  set-reset  latch,"  Appl.  Phys. 
Lett.,  52, 1419(1988). 

11.  M.E.  Prise,  N.  Striebl,  and  M.M.  Downs,  "Optical  considerations  in  the  design  of  a 
digital  optical  computer,"  Optical  and  Quantum  Electronics,  20, 49-77, 1988. 

12.  M.E.  Prise,  M.M.  Downs,  F.B.  McCormick,  S  J.  Walker,  and  N.  Streibl,  "Design  of  an 
optical  digital  computer,"  OSA  Topical  Meeting  on  Optical  Bistability,  OB-4  (March  1989). 

13.  M.E.  Prise,  N.C.  Craft,  R.E.  LaMarche,  M.M.  Downs,  S.J.  Walker,  L.A.  D’Asaro,  and 
L.M.F.  Chirovsky,  "Module  for  optical  logic  circuits  using  symmetric  self-electrooptic  effect 
devices,"  Appl.  Optics  29, 2164-2170  (1990). 

14.  A.  Dickenson  and  Michael  E.  Prise,  "Free-space  optical  intercoimection  scheme,"  Appl. 
Optics  29, 20012005  (1990). 

15.  M.  Banning,  "Practical  methods  of  making  and  using  multilayer  filters,"  J.  Opt.  Soc. 

Am.  37, 792-797  (1947). 

16.  J.  Mouchart,  J.  Begel  and  E.  Duda,  "Modified  MacNeille  cube  polarizer  for  a  wide 
an^lar  field,"  Appl.  Optics  28, 2847-2853  (1989). 

17.  R.A,  Chipman,  "Polarization  analysis  of  optical  systems,"  Opt.  Eng.  28,  90-999  (1989). 


73 


18.  R.A.  Chipman,  "Polarization  analysis  of  optical  systems  II,"  Proc.  SPIE  1166, 79-94 
(1989). 

19.  J.P.  McGuire,  Jr.  amd  R.A.  Chipman,  "Polarization  aberrations  in  optical  systems," 
Current  Developments  in  Optical  Engineering  II,  R.  Fischer,  W.Smith,  Editors,  Proc.  SPIE 
81, 240  (1987). 

20.  R.A.  Chipman  and  L.J.  Chipman,  "Polarization  aberration  diagrams,"  Opt.  Eng.  28, 
100-106  (1989). 

21.  J.L.  Pezzaniti  and  Russell  A.  Chipman,  "Imaging  Polarimeters  for  optical  system 
metrology,"  Proc.  SPIE  1317, 280-294  (1990). 

22.  M.F.  Morgan  and  R.A.  Chipman,  "An  ultraviolet  polarimeter  for  characterization  of  an 
imaging  spectrometer,"  Proc.  SPIE  1317,384-394  (1990). 


74 


Appendix  A 


Draft  of  a  paper  published  in  SPIE  proceedings  1317 


Imaging  polarimeters  for  optical  metrology 


J.  Larry  Pezzaniti 
Russell  A.  Chipman 


Dept,  of  Physics 

University  of  Alabama  in  Huntsville 
Huntsville,  AL  35899 


75 


Imaging  Polarimeters  for  Optical  Metrology 
Joseph  L.  Pezzaniti,  Russell  A.  Chipman 

University  of  Alabama  in  Huntsville,  Department  of  Physics 
Huntsville,  AL  35899 


Abstract 

In  this  paper  we  describe  two  configurations  of  imaging  polarimeters  designed  for 
polarimetric  optical  metrology.  The  first  is  a  Stokes  imaging  polarimeter  which  measures 
the  polarization  response  of  optical  systems  to  spherical  or  planar  waves  of  known 
polarization.  The  output  is  images  of  the  degree  of  polarization,  orientation  and 
eccentricity  of  polarization  ellipses,  or  Stokes  parameters  displayed  as  a  function  of  either 
the  exit  pupil  or  image  coordinate  of  the  optical  system.  The  second  configuration  is  a 
Mueller  imaging  polarimeter  which  measures  the  Mueller  matrix  of  an  optical  system  on  a 
ray  by  ray  basis.  Calibration  issues  involved  in  building  these  instruments  are  addressed 
along  with  a  brief  discussion  on  polarization  aberration  mechanisms. 


1.  Introduction 

For  the  majority  of  optical  systems,  detailed  information  about  their  effects  on  the 
polarization  of  light  is  not  important.  Systems  can  often  be  accurately  analyzed  by  ray 
tracing  using  the  assumption  that  polarization  and  intensity  across  the  wavefront  remains 
unaffected  by  the  optical  system.  When  this  assumption  is  no  longer  valid,  the  ray  tracing 
algorithm  must  be  generalized  to  include  the  polarization  effects  of  the  optical  system. 
Polarization  ray  tracing  calculates  the  variation  in  the  intensity  and  polarization  associated 
with  different  ray  paths  through  an  optical  system.i-'^  Similarly,  an  imaging  polarimeter  can 
be  used  as  an  optical  metrology  tool  to  measure  the  polarization  properties  of  optical 
systems  and  align  these  systems.  An  immediate  application  is  to  develop  a  data  base  on  the 
relative  properties  of  polarization  components  as  a  function  of  wavelength  and  field  of 
view. 

There  are  many  optical  systems  where  instrumental  polarizatiomis  a  concern,  such  as 
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interferometers,  polarimeters,  optical  computers,  grazing  incidence  optics,  and  phased 
arrays.  In  interferometric  systems,  such  as  laser  radars  and  analog  optical  computers,  light 
coupled  into  the  orthogonal  state  due  to  instrumental  polarization  will  not  interfere  with 
the  reference  beam,  decreasing  the  interferometer’s  fringe  visibility  and  information 
processing  capability.  Polarization  state  should  be  strictly  maintained  in  polarimeters  and 
ellipsometers,  since  any  spurious  polarization  contributions  from  the  instruments 
polarization  leads  to  inaccurate  results.  Other  examples  of  systems  in  which  polarization  is 
a  concern  include  broad  spectral  band  systems,  precision  radiometers,  phased  arrays,  and 
optical  signal  processors. 

The  University  of  Alabama  in  Huntsville  (UAH)  imaging  polarimeter  is  specifically 
designed  to  analyze  a  digital  optical  computer  under  development  at  AT&T  Bell 
Laboratories.5-7  This  system  uses  polarization  beam  combining  to  interlace  two  arrays  of 
beams  onto  a  modulator  and  bring  them  out  together  in  a  third  direction.  Each  logic 
module  uses  ten  polarizing  beam  splitters  and  twelve  quarter  wave  linear  retarders  and 
involves  fourteen  changes  of  polarization  state.  Small  polarization  effects  in  one  polarizing 
beam  splitter  can  become  serious  problems  in  a  series  of  ten  polarizing  beam  splitters. 
Furthermore,  systems  considerations  drive  the  designer  to  use  as  large  a  field  of  view  as 
possible  through  the  system.  If  the  polarization  state  is  not  correct,  then  at  each  polarizing 
beam  splitter,  part  of  the  beam  takes  the  wrong  path  and  eventually  is  imaged  onto  the 
wrong  modulator.  The  imaging  polarimeter  is  needed  for  the  study  and  alignment  of  this 
system. 

This  paper  describes  two  closely  related  types  of  imaging  polarimeters  for  metrology  of 
optical  systems,  a  Stokes  imaging  polarimeter  and  a  Mueller  imaging  polarimeter.  A 
Stokes  imaging  polarimeter  measures  the  polarization  response  of  an  optical  system  to  a 
given  input  polarization  state,  by  determining  the  Stokes  parameters  across  a  wavefront.  A 
Mueller  imaging  polarimeter  measures  the  Mueller  matrix  of  an  optical  system  on  a  ray  by 
ray  basis.  From  the  Mueller  matrix,  the  linear  and  circular  diattenuation  and  linear  and 
circular  retardance  associated  with  each  ray  path  may  be  calculated. 

2.  Stokes  imaging  polarimeter 

The  imaging  polarimeter  at  UAH  is  being  developed  for  optical  system  metrology. 
Other  Stokes  imaging  polarimeters  have  been  described  for  application  to  solar  vector 
magnetograph  measurements  8,  earth  resource  investigations  9,  and  remote  sensing.^O 

Figure  1  shows  a  schematic  of  our  Stokes  imaging  polarimeter.  It  consists  of  a  rotating 
retarder  polarimeter  followed  by  a  CCD  camera  coimected  to  an  image  processor.  The 
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Figure  1.  Schematic  diagram  of  a  rotating  retarder  Stokes  imaging  polarimeter. 
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Figure  2.  Schematic  diagram  of  a  Mueller  matrix  imaging  polarimeter. 
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polarimeter  consists  of  a  quarter  wave  linear  retarder  (QWLR)  followed  by  a  linear 
polarizer  (LP).  The  CCD  camera  has  a  576x485  resolution.  The  DT2861  image  processor 
board  holds  16  images  in  memory  simultaneously  and  is  able  to  do  frame 
addition/subtraction,  multiplication/  division,  zoom,  pan  and  scroll,  AND,  OR  ,XOR  logic 
operations,  and  display  output  in  false  color. 

The  objective  of  the  Stokes  imaging  polarimeter  is  to  determine  a  Stokes  vector 
associated  with  each  pixel  of  the  camera.  The  imaging  polarimeter  needs  to  acquire  a 
minimum  of  four  images  at  different  orientations  of  the  QWLR  to  calculate  the  Stokes 
parameters.  The  result,  following  data  reduction,  is  a  set  of  Stokes  images  where  a  Stokes 
image  is  defined  as  an  image  displaying  one  of  the  Stokes  parameters.  The  four  resulting 
Stokes  images  describe  the  intensity,  linear  polarization,  orientation,  and  circular 
polarization  content  of  the  light  in  the  scene. 

Our  interest  is  principally  in  polarimetric  metrology  of  optical  systems.  We  operate 
the  Stokes  imaging  polarimeter  by  introducing  a  plane  wave  or  spherical  wave  of  known 
polarization  state  into  the  optical  element  or  optical  system  under  test,  and  measuring 
Stokes  images  of  the  transmitted  light.  These  Stokes  images  can  be  measured  at  either  the 
exit  pupil  or  image  plane  by  focusing  the  camera  on  that  plane.  If  the  camera  is  focused 
on  the  exit  pupil,  then  each  pixel  of  the  camera  corresponds  to  a  ray  path  through  the 
optical  system  under  test  and  the  imaging  polarimeter  measures  the  polarization  aberration 
function.2  If  the  c^era  is  focused  on  the  image  plane,  then  the  distribution  of  polarization 
states  in  the  point  spread  function  is  measured  and  the  polarization  point  spread  function 
of  the  system  is  determined  in  Stokes  vector  form.  Since  the  point  spread  function  is 
usually  only  a  few  microns  in  diameter,  a  microscope  objective  or  other  high  magnific.  ion 
system  is  needed  to  observe  the  structure  of  the  point  spread  function  with  the  CCD 
camera. 

We  refer  to  pixel  coordinates  as  coordinates  which  locate  a  particular  pixel  on  the  CCD. 
These  pixel  coordinates  may  correspond  to  either  pupil  coordinates  or  image  coordinates 
of  the  optical  system  under  test,  depending  on  which  plane  the  imaging  polarimeter  is 
focussed. 

The  polarimeter  section  of  the  instrument  consists  of  a  QWLR  which  rotates  followed 
by  a  LP  which  is  fixed.  The  polarization  state  transmitted  through  the  polarimeter 
5 '  ( .  y ' )  *  written  as  a  function  of  pixel  coordinates  ( x ' ,  y ' ) ,  is  related  to  the 
polarization  state  incident  on  the  polarimeter  S(x ,  y ) ,  written  as  a  function  of  exit  pupil 
or  image  plane  coordinates  (x ,  y )  of  the  optical  system  under  test,  by  the  polarimeter 
equation 
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M  /./>(<!>)  is  the  Mueller  matrix  of  the  LP,  and  M  Qi/i;?  ( 0 )  is  the  Mueller  matrix  of  the  QWLR 
oriented  at  angle  0 .  The  LP  remains  at  a  fixed  orientation  (j)  =  0  for  all  measurements.  If 
the  linear  polarizer  changed  orientation,  the  instrumental  polarization  of  the  camera  lens 
and/or  CCD  might  need  to  be  compensated  as  a  possible  systematic  error.  This  would 
occur  if  the  transmission  through  the  lens  (or  a  folding  mirror)  or  the  sensitivity  of  the 
CCD  was  a  function  of  incident  polarization  state.  But  the  polarization  state  transmitted  to 
the  camera  lens  and  CCD  is  always  in  a  fixed  state  of  linear  polarization  insuring  that  the 
effect  of  the  instrumental  polarization  of  the  camera  lens  and  CCD  are  the  same  for  all 
measurements. 

For  the  LP  at  an  angle  (j)  =  0°  and  assuming  ideal  polarization  elements,  Eq.  1 
becomes: 
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The  CCD  measures  the  intensity  Sq  '  but  not  the  polarization  state  S'  of  tli.  'ieht. 

So '  is  a  function  of  the  QWLR  orientation  0, 

So'(.v'.y'.0)  =  i(2so(.Y.y)  +  s,(.v.y)(  1  +  cos40)  +  SjC.v.  y)sin  40- ZsjC.v.  y )sin 20 

(3) 

Video  frames  of  Sq  ' (.v ' .  y ' ,  0)  are  digitized  at  four  different  orientations  0  of  the 
QWLR.  By  orienting  the  QWLR  at  orientations  that  produce  four  linearly  independent 
equations  in  So(.v,y),Si(x .  y).S2(,v,y),S3(x.y),  the  four  Stokes  images  are 
calculated  from  the  four  measured  images  So'(.v.y,0,),So'(.v,y,02), 
So'(x,y,Qj),SQ'ix,y,Q^). 

The  QWLR  is  necessary  to  obtain  the  third  Stokes  parameter  S3  because  a  polarimeter 
consisting  of  only  a  LP  can  measure  only  Sq  ,  s , ,  Sa . 
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One  data  acquisition  routine  which  provides  a  particularly  simple  set  of  inversion 
equations  orients  the  QWLR  at  angles  0  =  45° .  -45° .  67.5° ,  90°  The  four  Stokes 
images  are  obtained  from  the  equations 


^o(^.y)  =  So' y'. -15°  )  +  -45°),  (4) 

s,(.v.y)  =  2so'(.v'.y'.90°)-So'(.v'.y',4S°)-5o'(.v'.y'.-45°).  (5) 

S2(.v.y)  =  (l+'/2)So'(-v'-y'.-15°)-(l-v2)So'(.v'.y'.-45°) 

■^2so'(.v'.y'.90°)-4So'(-^-'.y'.67.5°).  (6) 

SaC'V.y)  =  So'(.v',y'.-45°)-So'(.v',y',45°).  (7) 


While  images  acquired  at  QWLR  angles  of  0  =  45° ,  -45° ,  67.5° ,  90° ,  simplify  the 
data  reduction,  they  do  not  measure  all  four  Stokes  parameters  with  equal  accuracy.  The 
first  two  measurements  at  0  ®=  45°  and  0  =  -  45°  involve  only  So  and  S3 ,  while  the 
measurement  at  0  =  67.5°  involves  so.s,  .Sj  and  S3,  and  0  =  90°  involves  only  soandsi . 
"Therefore  So  and  S3  are  determined  from  two  measurements,  s  1  is  determined  from  three 
measurements,  and  sg  is  determined  from  four. 

Two  different  methods  have  been  developed  to  calculate  the  Stokes  images.  In  the  first 
method  the  So  '(0 .  .v ,  y )  frames  are  added  and  subtracted  in  the  Data  Translation 
DT2861  frame  grabber  board,  according  to  Eqs.  4-7.  The  accuracy  of  this  method  is 
restricted  by  the  DT2861  frame  grabber  board  since  it  performs  integer  arithmetic  with  8 
bits.  If  an  addition,  multiplication  or  subtraction  results  in  a  pixel  level  above  255  or  below 
0,  the  value  will  become  0  or  255.  For  this  reason  levels  in  the  Sq  ' ( 0 ,  -v ,  y )  image 
contrast  must  be  compressed  by  a  factor  of  at  least  four  with  a  corresponding  loss  of 
accuracy. 

The  second  method  performs  the  Stokes  vector  calculations  as  floating  point 
calculations  in  computer  memory.  'This  does  not  require  contrast  compression  but  takes 
longer  to  calculate.  The  arrays  containing  the  Stokes  images  are  then  shipped  back  over 
the  computer  bus  to  the  frame  buffers  of  the  DT2861  board,  where  the  images  are 
displayed.  Presently,  the  four  Stokes  vector  images  are  calculated  in  about  15  sec,  by 
method  1  and  about  3  min.  by  method  2. 

The  Stokes  images  contain  intensity  as  well  as  polarization  information,  all  as  a  function 
of  image  coordinate.  The  zeroth  Stokes  image  So  y )  is  the  intensity  image  and  is 
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always  positive.  The  first,  second  and  third  Stokes  images  can  be  positive  or  negative,  so 
the  images  are  displayed  with  zero  as  128.  Positive  values  of  s  x,  s  2,  are  displayed  from 
129  to  255,  and  negative  values  from  0  to  127,  requiring  a  contrast  compression  of  two. 

A  number  of  other  useful  measures  of  polarization  can  be  displayed  as  images.  The 
degree  of  polarization  DOP  is  a  measure  of  the  overall  polarization  of  the  light.  It  is 
defined  as 


DOP  = 


.  ->  2 


o2 

3 


DOP  varies  from  zero  for  unpolarized  light  to  one  for  polarized  light. 

The  degree  of  linear  polarization  DOLP  is  a  measure  of  the  amount  of  linear 
polarization 


(8) 


DOLP  = 


JO 


(9) 


DOLP  is  zero  for  unpolarized  or  circularly  polarized  light  and  is  one  for  linearly  polarized 
light. 

The  degree  of  circular  polarization  DOCP  is  a  measure  of  the  amount  of  circular 
polarization 


IS3 1 

DOCP  =  ^-^. 
So 


(10) 


DOCP  is  zero  for  unpolarized  and  linearly  polarized  light  and  one  for  right  or  left  circularly 
polarized  light. 

Two  other  definitions  are  necessary  to  specify  a  polarization  ellipse.  The  eccentricity  x 
of  the  polarization  ellipse  is 


X  =  -arcsin\ 


si 


VsT  +  s|*s5. 


and  the  orientation  a  of  the  major  axis  of  the  polarization  ellipse  is 


(11) 
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a  =  ^arctani—\  (12) 

2  \sj 

One  useful  display  of  polarization  aberration  is  graphical  maps  containing  arrays  of 
polarization  ellipses  representing  the  average  polarization  of  regions  of  the  wavefront,  such 
as  is  shown  in  Figure  3. 

Similarly  images  of  DOP(x ,  y) ,  DOLP(x .  y  )and  DOCP(.v ,  y  )can  be  displayed  with 
256  grey  levels  or  in  false  color. 

3.  Mueller  matrix  imaging  polarimeter 

Mueller  matrix  polarimeters  have  been  built  to  measure  the  Mueller  matrix  of 
polarizing  elements  as  a  function  of  wavelength>li  to  measure  Mueller  matrices  at  laser 
wavelengths, 12  and  to  measure  the  Mueller  matrices  of  scattered  light.13,14  The  UAH 
Mueller  matrbc  imaging  polarimeter  will  measure  the  Mueller  matrix  as  a  function  of  exit 
pupil  coordinate  or  image  coordinate  of  an  optical  system. 

The  Mueller  matrices  describing  optical  systems  are  in  general  functions  of  object  and 
pupil  coordinate.  Each  ray  path  through  an  optical  system  will  have  slightly  different 
polarizing  properties  because  of  different  angles  of  incidence  and  orientations  of  its  plane 
of  incidence  at  thin  films  and  different  directions  and  path  lengths  through  birefringent  or 
dichfoic  materials.  The  polarization  aberration  function  iVl(H.p,\)  for  an  optical  system 
is  a  Mueller  matrbc  as  a  function  of  object  coordinate  h.,  pupil  coordinate  p  ,  and 
wavelength  A.  1.2  The  spatial  dependence  of  the  Mueller  matrix  gives  rise  to  polarization 
aberrations  at  the  exit  pupil  of  the  optical  system.  1-4  The  Mueller  matrbc  imaging 
polarimeter  measures  the  spatial  dependence  of  the  Mueller  matrbc  for  a  set  of  ray  paths 
through  an  optical  system  under  test. 

The  UAH  imaging  polarimeter  uses  a  collimated  quasimonochromatic  light  source  and 
two  QWLR  -  LP  pairs  in  conjunction  with  a  CCD  camera  (see  figure  2).  Each  pbcel  of  the 
CCD  camera  measures  the  Mueller  matrbc  for  a  particular  ray  path  through  the  optical 
system.  The  light  from  the  source  is  collimated  and  passes  through  a  LP  at  a  fbced 
orientation  angle,  then  a  QWLR  mounted  to  a  rotary  stage.  The  light  source  and  the 
LP/QWLR  pair  comprise  the  polarization  generator.  Lenses  then  prepare  the  appropriate 
input  spherical  wave  for  the  optical  system  under  test,  where  the  polarization  of  each  ray 
will  be  altered  according  to  the  Mueller  matrbc  for  that  particular  ray  path.  The  light  is 
recollimated  and  enters  the  Stokes  imaging  polarimeter.  The  camera  is  focused  either  on 
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the  exit  pupil  of  the  optical  system  or  its  image  plane.  By  fixing  the  orientation  of  the  two 
linear  polarizers,  the  imaging  polarimeter  becomes  insulated  from  instrumental 
polarization  before  the  first  polarizer  and  after  the  final  polarizer. 

Any  optics  between  the  polarizers  used  to  reshape  the  wavefront  may  introduce 
additional  instrumental  polarization.  This  potential  systematic  error  must  be  compensated 
for  in  the  calibration  and  data  reduction  if  the  polarization  contribution  from  these 
elements  are  substantial. 

Our  Mueller  matrix  imaging  polarimeter  will  use  a  variation  of  the  method  proposed  by 
.4zzam  15  and  expanded  by  others  16,17  to  calculate  the  Mueller  matrices  from  the  raw 
polarimetric  data.  In  Azzam’s  method,  the  orientation  of  the  two  linear  polarizers  are  fixed 
and  a  series  of  at  least  24  images  are  recorded  with  different  retarder  settings.  The  first 
image  is  acquired  with  the  transmission  axes  of  the  polarizers  and  the  fast  axes  of  the 
retarders  all  parallel.  Between  image  measurements,  the  first  retarder  is  rotated  by  A0 
and  the  second  retarder  by  5  A0  until  the  first  retarder  has  rotated  through  180°.  For 
example,  if  the  generating  retarder  is  rotated  in  increments  of  5  degrees,  the  analyzing 
retarder  is  rotated  in  25  degree  increments.  This  generates  an  intensity  modulation  at  each 
pixel  in  the  images.  The  polarization  information  is  encoded  on  the  first  twelve  harmonics 
and  the  dc  component,  so  a  minimum  of  twenty  four  images  are  needed  over  180  degrees 
to  sample  the  twelfth  harmonic.  The  Fourier  coefficients  are  calculated  with  the  discrete 
Fourier  transform  at  each  pbcel.  The  Fourier  coefficients  are  manipulated  to  obtain  the 
sixteen  elements  of  the  Mueller  matrix  associated  with  each  pixel. 

The  Mueller  matrix  can  be  displayed  as  images  of  the  individual  matrix  elements,  or  in 
other  related  polarization  output.  Our  principal  interest  is  to  display  the  linear  retardance 
and  linear  diattenuation  in  image  form. 

4.  Calibration  issues  of  an  imaging  polarimeter 

Imaging  polarimeters  which  take  accurate  data  are  very  difficult  to  build.  It  is  easy  to 
buy  the  components  and  assemble  an  instrument.  It  is  a  long  arduous  project  to  align  and 
calibrate  the  polarimeter,  write  software  to  acquire  and  reduce  the  data,  and  address 
meaningful  optical  metrology  problems.  In.this  section,  the  issues  of  bootstrap  calibration, 
drift,  moving  objects,  element  wedge,  angle  of  inciaence  dependence  in  the  polarization 
elements,  and  instrumental  polarization,  are  discussed. 

If  the  Mueller  matrices  for  the  polarizers  and  retarders  used  in  the  polarimeter  differ 
from  their  ideal  values  in  Eq.  2,  then  the  data  reduction  routines  (Eq.  4-7)  should  be 
modified  to  take  their  actual  behavior  into  account.  Commercially  available  polarizers  and 
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retarders  are  not  specified  well  enough  by  the  vendors  to  accurately  estimate  their  Mueller 
matrices,  nor  are  instruments  or  test  services  readily  available  to  measure  these  Mueller 
matrices.  Retarders  are  never  exactly  a  quarter  wave.  Polarizers  always  have  some 
leakage.  Retarder  and  polarizer  characteristics  vary  with  wavelength  and  angle  of 
incidence.  Retarders  display  some  diattenuation.  Polarizers  display  some  retardance. 
Linear  polarizers  or  linear  retarders  may  actually  be  elliptical. 

Without  sufficient  prior  knowledge  of  the  properties  of  the  polarization  elements,  a 
bootstrap  calibration  needs  to  performed,  i.  e.  a  calibration  started  using  retarders  and 
polarizers  with  unknown  properties.  The  calibration  should  make  the  fewest  number  of 
assumptions  about  the  polarizing  elements  and  utilize  methods  that  are  relatively 
insensitive  to  the  weakest  assumptions.  Once  an  accurate  Mueller  matrix  polarimeter  is 
operating,  it  can  be  used  to  calibrate  polarization  elements  for  other  polarimeters,  greatly 
simplifying  the  exacting  task  of  construction  an  accurate  imaging  polarimeter. 

Preliminary  calibration  has  been  done  on  the  UAH  imaging  Stokes  polarimeter. 
Polarizer  calibration  assumed  that  the  polarizers  only  display  diattenuation.  The 
transmittance  of  the  transmission  axis  and  the  absorption  axis  of  polarizers  were  measured. 
Two  similar  polarizers  were  calibrated  in  this  way  were  then  used  to  calibrate  the  retarders. 
The  first  nine  elements  rn  1 1  through  rrij^  of  the  retarder  Mueller  matrix  were  measured 
by  rotating  the  polarizers  at  a  3:1  ratio  on  either  side  of  the  QWLR.  The  resulting  intensity 
modulation  has  four  harmonics  and  a  dc  term  in  its  Fourier  series  which  will  solve  for  the 
nine  Mueller  matrix  elements.  The  linear  retardance  and  linear  diattenuation  can  be 
obtained  from  these  nine  Mueller  matrbe  components  assuming  there  is  no  circular 
retardance.  Circular  retardance  does  occur  in  multiple  element  linear  retarders  when  the 
fast  axes  of  the  linearly  birefringent  elements  are  not  exactly  0°  or  90°  apart.  To  test  for 
circular  retardance,  the  QWLR  was  placed  between  two  crossed  linear  polarizers  with  its 
fast  parallel  to  the  transmission  axis  of  the  first  polarizer.  If  circular  retardance  is  present, 
the  electric  field  vector  is  rotated  slightly,  so  that  the  final  polarizer  must  be  rotated  to 
locate  the  null.  No  circular  retardance  was  observed  in  these  retarders,  but  we  have 
detected  it  in  other  linear  retarders.  If  circular  retardance  is  observed,  the  remaining  seven 
components  of  the  QWLR  Mueller  matrbe  should  be  measured.  A  method  to  measure 
these  components  is  described  in  reference  18. 

The  next  issue  is  drift.  The  imaging  polarimeters  discussed  here  take  measurements 
sequentially  in  time.  Any  change  in  source  intensity  or  any  drift  in  camera  sensitivity 
during  the  measurement  sequence  is  interpreted  as  polarization,  albeit  spurious 
polarization. 

Similar  problems  arise  in  polarization  measurements  of  a  moving  object.  Consider  how 
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digital  image  processors  calculate  derivatives.  A  partial  derivative  of  the  image  intensity  in 
the  X  direction  is  calculated  by  shifting  the  image  in  the  x  direction  and  subtracted  the 
image  from  itself.  The  imaging  polarimeter  performs  data  reduction  by  adding  and 
subtracting  images.  A  moving  object  in  the  scene  introduces  spurious  contributions  into 
the  Stokes  images  proportional  to  the  resulting  intensity  derivative.^.  These  intensity 
gradients  are  often  greatest  around  the  objects  edges,  but  edges  are  where  many  objects 
have  their  greatest  polarization  signatures.  Fortunately,  in  optical  system  polarimetric 
metrology,  our  subject  is  usually  a  fixed  optical  system  and  the  intensity  gradients  are 
generally  small.  Thus,  far  more  accurate  polarimetry  is  possible  in  this  setting  than  in 
remote  sensing  or  tactical  polarimetry. 

The  faces  of  the  QWLR  need  to  be  highly  parallel  because  wedge  in  a  rotating  element 
generates  beam  wander.  If  the  retarder  is  slightly  prismatic,  the  image  will  move  in  a  circle 
on  the  CCD  resulting  in  a  spatial  misregistration  and  image  derivatives.  Large  amounts  of 
beam  wander  can  be  partially  compensated  by  carefully  shifting  the  measured  images  in 
software,  but  this  correction  cannot  be  performed  at  the  fractional  pixel  level. 

If  light  propagates  at  substantial  angles  through  the  polarization  generator  and  the 
polarization  analyzer,  the  field  angle  dependence  of  the  polarizers  and  the  retarders  has  to 
be  treated  as  systematic  error.  For  exit  pupil  measurements,  each  pixel  on  the  CCD  will 
corresponds  to  a  different  angle  through  the  polarization  components.  A  straightforv»'ard 
compensation  can  be  performed  which  used  the  Mueller  matrices  as  a  function  of 
propagation  direction  in  the  data  reduction. 

Tne  instrumental  polarization  due  to  wavefront  shaping  optics  within  the  polarimeter 
may  require  compensation  during  data  reduction.  Their  instrumental  polarization  will 
have  to  be  determined  each  time  a  new  optical  system  is  configured  within  the  sample 
compartment,  so  standardized  methods  for  removing  this  systematic  error  should  be 
developed. 


5.  lypical  polarization  aberration  mechanisms 

The  objective  of  the  UAH  imaging  polarimeter  is  to  study  polarization  aberration  in 
optical  components  and  optical  systems.  Polarization  aberration  contributions  from  optical 
elements  can  be  conveniently  grouped  into  two  categories:  instrumental  polarization  and 
nonideal  polarization.  The  first  category,  "instrumental  polarization"  is  due  to  optical 
elements  which  are  not  intended  to  control  polarization  such  as  lenses,  mirrors  and  prisms. 
These  elements  usually  have  weak  polarization.  The  other  category,  "nonideal 
polarization,"  is  the  deviation  in  performance  of  elements  used  to  control  polarization 
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from  itteir  ideal  characteristics.  Several  examples  of  polarization  aberration  problems  are 
mentioned'here,  microscope  objectives,  mirrors,  coating  thickness  variations,  field 
dependence  in  retarders,  and  polarizing  beam  splitters. 

The  classic  example  of  polarization  aberration  is  the  instrumental  polarization  of 
microscope  objectives  in  polarization  microscopy.  20  jhe  Fresnel  transmission 
coefficients  cause  a  lens  to  act  with  a  radially  oriented  linear  diattenuation  which  increases 
quadratically  from  the  center  of  the  pupil.  This  causes  a  rotation  of  the  plane  of 
polarization  in  four  corners  of  the  pupil  when  linearly  polarized  light  is  incident.  The 
objective  lens  of  a  microscope  contributes  a  substantial  leakage  in  a  "Maltese  cross"  pattern 
when  located  between  crossed  polarizers. 

Metal  mirrors  introduce  linear  retardance  due  to  the  difference  between  s  and  p  phase 
shifts  upon  reflection  along  with  a  lesser  amount  of  linear  diattenuation.  Figure  3  shows  a 
polarization  ray  trace  calculation  of  a  f/0.5  parabolic  mirror  coated  with  aluminum 
illuminated  with  circularly  polarized  incident  light.  At  the  center  of  the  mirror,  where  the 
rays  are  normally  incident,  there  is  no  difference  in  phase  shift  between  s  and  p  light.  The 
handedness  of  the  circular  polarization  changes  upon  reflection.  At  the  edge  of  the  mirror 
the  linear  retardance.approaches  a  quarter  wave  with  the  effect  of  coupling  circular 
polarized  light  to  nearly  linearly  polarized  light  oriented  at  45 “to  the  p  plane.  Another 
example  of  instrumental  polarization  from  mirrors  is  the  very  large  polarization  effects  \vith 
grazing  incidence  telescopes  such  as  the  Advanced  X*ray  Astrophysical  Facility  (AXAF) 
telescope  under  development  by  NASA.21 

Polarization  aberration  results  from  variations  in  multi-layer  thin  film  coatings 
thicknesses.  For  instance,  with  coatings  of  ten  to  thirty  layers,  the  retardance  can  vary  tens 
c.f  degrees  over  an  f/2  parabolic  mirror  due  to  coating  thickness  variations  and  angle  of 

incidence  variations.22 

The  second  category  of  polarization  aberration  is  nonideal  polarization  from  elements 
designed  to  control  polarization  including  polarizers,  retarders  and  polarizing  beam 
splitters. 

Several  characteristics  of  birefringent  retarders  introduce  polarization  aberration  in 
optical  systems.  First,  the  thickness  of  the  birefringent  material  may  be  wrong  or  even 
wedged,  yielding  an  incorrect  or  spatially  varying  retardance.  The  crystal  axis  can  be 
misoriented  relative  to  the  faces  or  the  mount.  The  retarder  can  be  misaligned  in  the 
optical  system,  either  tilted  or  with  its  fast  axis  at  the  wrong  orientation. 

Three  other  problems  with  retarders  are  more  fundamental,  field  dependence, 
diattenuation,  ..and  birefraction. 

Figure  4  shows  how  the  extraordinary  refractive  index  and  the  retardance  varies  with  the 
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direction  of  the  incident  light.  Different  parts  of  a  spherical  wave  propagating  through  a 
bireffingent  retarder  experience  different  retardances  yielding  polarization  aberr:  t  .on. 

Retarders  have  some  weak  diattenuation,  both  in  their  bulk  absorption  spectrua:  and 
from  the  difference  in  Fresnel  losses  between  the  ordinary  and  extraordinary  rays  at  the 
interfaces.  Anti-reflection  coatings,  carefully  selected,  (an  greatly  reduce  the  diattenuation 
from  birefringent  interfaces. 

Figure  5  shows  the  birefraction  or  double  refraction  in  birefringent  retarders  and 
uniaxial  media.  The  crystal  axis  is  represented  by  c.  If  the  propagation  vector  in  the 
crystal  is  not  parallel  or  perpendicular  to  the  crystal  axis  of  the  uniaxial  media,  the  ordinary 
and  extraordinary  rays  take  divergent  paths.  This  causes  two  crescents  on  the  edges  of  the 
beam  to  have  different  polarization  states  from  the  main  beam.  Fig.  5  shows  two  cases  of 
illumination  of  a  quarter  wave  retarder  at  different  angles  of  incidence  with  beams  linearly 
polarized  at  45**  with  respect  to  the  fast  axis.  In  the  first  case,  the  beam  is  normally 
inddent  and  ^e  entire  transmitted  beam  is  circularly  polarized.  In  the  bottom  example, 
the  beam  is  not  normally  inddent.  Since  the  refractive  indices  for  the  extraordinary  and 
ordinary  rays,  represented  by  the  lines  and  dots  respectively,  are  different,  the  two  rays  are 
refracted  at  different  angles.  The  transmitted  beam  has  the  desired  polarization  where  the 
two  beams  overlap.  The  aescent-shaped  areas  along  either  edge  of  the  light  bundle 
contain  light  polarized  in  the  eigen-polarization  states  of  the  medium. 

Polarizing  beam  splitters  are  notorious  for  causing  systems  problems  due  to  nonideal 
polarization.  A  ideal  polarizing  beam  splitter  reflects  all  the  s  polarizeci  light  and  transmits 
all  the  p  polarized  light.  In  many  applications,  such  as  digital  optical  computers,  a  freld  of 
view  of  many  degrees  is  required,  so  the  reflection  and  transmission  coefficients  should  be 
nearly  constant  over  this  range  of  field  angles.  Since  the  light  is  at  non-normal  incidence  at 
the  beam  splitting  interface,  the  reflection  and  transmission  coefficients  can  vary  linearly  in 
the  angle  of  incidence.  This  makes  polarizing  beam  splitters  intrinsically  more  angle  of 
incidence  sensitive  than  mirror  or  lens  coatings  which  can  show  variation  to  only  second 
order  about  normal  incidence.  We  have  yet  to  see  a  good  polarizing  beam  splitter  and  are 
interested  if  there  are  fundamental  limitations  on  the  field  of  view  performance  of 
polarizing  beam  splitters. 


12 


Flgarc  5.  figure  shows  how  beam  shear  between  extraordinary  and  ordicary  rays  lim-its  the  held  of  view  of 
a  quarter  wa  ve  retarder. 


Imaging  Polarimeters  for  Optical  Metrology 


6.  Conclusion 

The  imaging  polarimeter  is  a  tool  to  measure  the  instrumental  polarization  of  optical 
systems  and  calibrate  polarization  elements.  It  measures  polarization  properties  associated 
with  ray  paths  through  optical  systems.  It  can  also  determine  the  polarization  state 
structure  of  an  image.  It  provides,  via  CCD  camera  and  image  processing  system,  graphical 
maps  of  polarization  information.  It  can  be  used  to  assist  in  the  alignment  of  polarization 
critical  optical  systems.  Finally  it  is  needed  for  comparative  data  on  the  performance  of 
polarizing  elements,  where  it  can  simultaneously  display  performance  as  a  function  of  field 
of  view  in  an  image  format. 

It  is  easy  in  principle  to  assemble  an  imaging  polarimeter,  but  quite  difficult  to  obtain 
accurate  results.  Each  of  the  calibr.ition  issues  is  time  consuming,  and  often  maddeningly 
coupled  into  each  other. 

Currently  the  UAH  imaging  polarimeter  is  able  to  measure  Stokes  images,  degree  of 
polarization  images,  degree  of  linear  'lOlarization  images,  degree  of  circular  polarization 
images,  and  images  showing  the  orientation  of  the  major  axis  of  a  polarization  ellipse.  A 
preliminary  calibration  has  been  done  on  the  QWLR  and  LP  at  633  nm.  Soon,  the 
polarimeter  will  measure  Mueller  matrices  at  633  nm  and  849  nm.  A  study  on  the 
performance  of  polarizing  beam  splitter  cubes  as  a  function  of  field  angle  is  in  progress. 

We  also  intend  to  study  depolarization,  the  coupling  of  polarized  light  into  unpolarized 
light,  in  optical  systems. 
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Angular  dependence  of  polarizing  beam  splitter  cubes 


Abstract 

A  method  to  measure  the  angular  properties  of  polarizing  beam  splitters  was  developed  for 
characterizing  polarizing  beam  splitters  used  in  polarization  based  imaging  systems. 

Several  performance  characteristics,  such  as  transmission,  reflection,  leakage,  and  variation 
of  the  orientation  of  linear  polarization  were  measured  with  an  imaging  polarimeter. 
.Vleasurements  were  made  over  a  10°  range  of  angle  of  incidence  simultaneously  on  three 
commercially  available  polarizing  beam  splitters,  and  one  specifically  designed  for  a  wide 
angular  field  of  view.  Variation  of  the  orientation  of  linear  polarization  of  polarizing  beam 
Npliiters  is  described.  Derivation  of  the  Mueller  matrbc  of  a  polarizing  beam  splitter  is 
provided  in  an  Appendbt. 

I.  Introduction 

An  ideal  polarizing  beam  splitter  divides  the  incident  light  into  s  and  p  polarization 
states,  transmitting  the  p  polarized  component  and  reflecting  the  s  polarized  component. 
Polarizing  beam  splitters  have  a  limited  range  of  angle  of  incidence  and  wavelength  for 
which  they  are  effective.  Typically,  polarizing  beam  splitter  designs  are  optimized  to 
provide  a  large  spectral  range  within  a  small  angular  field  of  view.  For  example,  the 
.Vlacneille  polarizing  beam  splitter  design  achieves  extinction  ratios  of  greater  than  100:1, 


biit.the  beam  must  be  collimated  to  within  ±2°  of  the  entrance  face  noxma\..[l] [2] 
Applications  which  involve  imaging  or  non-collimated  light,  however,  may  require  the 
polarizing  beam  splitter  to  have  a  wide  field  of  view,  often  with  quasimonochromatic  light. 

The  angular  behavior  of  polarizing  beamsplitters  is  the  subject  of  this  paper.  The 
objective  is  to  characterize  polarizing  beam  splitter  behavior  in  polarization  based  imaging 
.systems,  specifically  optical  computers.  For  example,  several  free  space  optical 
interconnects  utilize  the  polarizing  beam  splitter’s  four  port  access  to  losslessly  interlace 
and  divide  arrays  of  beams  carrying  logic  \TdoTm2L\xovi.[3] [4] [5] [6] [7] [8]  In  these  systems, 
2-D  logic  arrays  are  interconnected  by  imaging  arrays  of  spots  (generated  by  binary  phase 
gratings)  from  one  device  to  the  next.  The  logic  devices  are  situated  in  image  planes  and 
the  polarizing  beam  splitters  are  located  in  planes  at  or  near  the  pupil  where  the  light  is 
collimated.  Collimated  light  from  each  element  of  a  logic  array  passes  through  the 
polarizing  beam  splitters  at  a  different  angle.  Presently,  the  above  referenced  systems 
o|)crate  at  =.T  I'iokl  of  view,  but  future  systems  will  be  extended  to  the  largest  fields  of  view 
permi.ssible.i  10'’  or  more.  Since  most  system  designs  require  each  beam  to  pass  through 
several  polarizing  beam  splitters,  the  system  efficiency  depends  on  the  angular  response  of 
the  polarizing  beam  splitters.  Furthermore,  light  incorrectly  divided  at  the  polarizing  beam 
splitter  will  usually  be  imaged  onto  an  incorrect  logic  element  or  detector. 

This  paper  introduces  phenomenological  data  taken  on  three  commercially  available 
polarizing  beam  splitters  and  one  polarizing  beam  splitter  specifically  designed  for  a  field 
of  view  of  Six  performance  parameters  are  measured  for  each  polarizing  beam 

splitter  using  an  imaging  polarimeter/iO/.  The  polarizing  beam  splitters  are  measured  in 
transmission  and  reflection  over  a  ±10°  field  of  view. 

One  aspect  of  polarizing  beam  splitters  examined  by  this  paper  is  the  dependence  of 
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orientation  of  linear  polarization  on  ray  angle  of  incidence.  In  Section  II  we  provide  a 
conci.se  explanation  of  how  and  why  the  orientation  of  linear  polarization  in  the  reflected 
and  transmitted  beams  depend  on  ray  angle  of  incidence.  Section  III  introduces  the  six 
figures  of  merit  which  we  measured  on  each  of  the  four  polarizing  beam  splitters.  Section 
IV  describes  our  imaging  polarimeter  and  the  procedures  used  to  measure  the  angular 
dependence  of  the  polarizing  beam  splitters.  Section  V  presents  the  data  on  the  polarizing 
beam  splitters.  Finally,  Appendbc  A  derives  the  Mueller  matrix  of  a  polarizing  beam 
.splitter  in  transmi.ssion  and  reflection  in  terms  of  the  complex  transmission  and  reflection 
coefficients  of  the  beamsplitting  interface,  relating  our  measurements,  the  Mueller 
calculus,  and  the  quantities  calculated  by  thin  film  analysis  codes,  s  and  p  reflectance, 
transmittance,  and  pha.se  difference. 

n.  Variation  of  the  Orientation  of  Linear  Polarization 

A  polarizing  beam  splitter  cube  is  formed  by  coating  the  hypotenuse  face  of  an 
iso.sceles  right-angled  prism  with  a  number  of  thin-film  layers  and  then  cementing  it  to  an 
identical  uncoated  prism.  The  usual  design  of  the  multi-layer  has  been  of  the  form  initially 
devised  by  Macneal  [1  ]  and  consists  of  alternate  layers  of  a  high  and  a  low  refractive  index 
material.  The  optical  thickness  of  each  layer  is  equal  to  a  quarter  of  the  design  wavelength 
at  the  angle  of  refraction  in  that  layer,  and  the  refractive  indices  of  the  thin-film  materials 
as  well  as  that  of  the  substrate  are  chosen  so  that  the  light  is  incident  upon  each  of  the 
multi-layer  at  the  Brewster  angle  [11] [12] [13].  Due  to  Brewster’s  condition  the 
reflectance  of  the  p  polarization  approaches  zero  while  the  reflectance  of  the  s  polarization 
remains  high.  Angles  of  incidence  different  from  45°,  however,  changes  the  optical  path 
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length  from  the  quarter  wave  configuration,  and  Brewster’s  condition  no  longer  holds. 

This  causes  a  leakage  of  s  polarization  in  transmission  and  p  polarization  in  reflection. 
References  f  141(15]  contain  two  different  approaches  to  polarizing  beam  splitter 
multi-layer  design  to  increase  the  angular  range  for  which  the  p  transmittance  and  s 
reflectance  remain  high. 

The  orientation  of  linear  polarization  transmitted  from  a  polarizing  beam  splitter  is 
dependent  on  the  ray  angle  of  incidence.  If  a  spherical  wave  of  unpolarized  light  is  incident 
on  a  polarizing  beam  splitter,  the  reflected  and  transmitted  beams  will  have  a  spatially 
varying  linear  polarization.  This  is  because  the  orientation  of  linear  polarization  of  a 
polarizing  beam  splitter  aligns  with  the  plane  of  incidence  (p  planej  in  transmission  and  the 
s  direction  in  reflection,  while  the  s  and  p  orientations  depend  on  the  ray’s  angle  of 
incidence. 

Figure  1  introduces  a  coordinate  system  for  a  polarizing  beam  splitter  cube,  '1  he 
origin  is  centered  on  the  entrance  face  with  the  z  axis  normal  to  the  face.  The  beam 
splitting  interface  normal  lies  in  the  x-z  plane  and  is  45°  from  the  z  axis.  The  direction  of 
an  incident  ray  is  specified  by  its  refracted  direction  cosines  (l.m.ri)  where  I  =  cosa , 

/?i  =  c{)o|5  and /t  =  cosy  and  i^  +  m^  +  n^=  1.  a, 3,  and  y  are  illustrated  in  Figure  1.  For 
rays  incident  such  that  3  =  90° ,  the  p  direction  is  parallel  to  the  x-z  plane  (we  will  refer  to 
this  as  the  horizontal  plane)  and  the  s  direction  aligns  with  the  y-z  plane  (vertical).  For 
angles  of  incidence  in  which  3  >  90°,  the  p  plane  rotates  counterclockwise  with  the  respect 
to  the  horizontal,  while  for  3  <  90°  it  rotates  clockwise.  Figure  2  illustrates  the  orientation 
of  s  and  p  directions  as  functions  of  the  direction  cosines  of  the  incident  rays.  In  Figure  2, 

■  increases  from  90°  in  the  positive  x  direction  and  3  increases  from  90°  in  the  positive  y 
tlirection.  The  thin  cross  hairs  are  parallel  with  the  horizontal  and  vertical  planes  (,x-z  and 
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y-z  planes  in  Figure  1).  The  large  circles  indicate  the  zones  where  the  angles  of  incidence 
are  y  =  5°  and  y=  10°.  The  angle  4)  is  the  angle  between  the  p  direction  and  the 
horizontal  x-z  plane. 

Consider  an  ideal  polarizing  beam  splitter  which  transmits  100%  of  the  p  polarized 
light  and  reflects  100%  of  the  s  polarized  light  for  all  angles  of  incidence.  Two  examples 
readily  illustrate  the  polarization  aberration  which  polarizing  beam  splitters  can  introduce 
in  imaging  systems:  1.  a  spherical  wave  of  linearly  polarized  light,  and  2.  a  polarizing  beam 
splitter  at  the  pupil  of  an  4F  optical  filtering  system.  For  the  first  case,  consider  light 
uniformly  polarized  parallel  to  the  x-z  plane  focused  through  a  polarizing  beam  splitter. 

The  situation  is  illustrated  in  Figure  3,  The  s-p  coordinates  rotate  for  rays  in  which 
p;  *  90° .  The  intensity  of  each  ray  transmitted  depends  on  the  angle  (j)  between  the  linear 
polarization  direction  of  the  incident  ray  and  the  p  direction  at  the  beamsplitting  interface. 
,'\ccorciing  to  Malus’slaw  the  transmitted  intensity  of  a  ray  is 

/(d))- /oCoij^'4)  (1) 

where  <P  Ls  a  function  of  the.incident  rays  direction  cosines,  and  /ois  the  incident  intensity. 
Eqn.  A I  i  in  Appendix  A  contains  an  analytic:expr.p.ssion  for  An  ideal 

polarizing  beam  splitter  cube  introduces  apodization  and  a  spatial  variation  of  linear 
polarization  orientation,  aligned  with  p  or  s  directions,  inuhe  transmitted  and  reflected 
beam, 

For  the  second  example,  consider  a  polarizing  beam  splitter  situated  between  infinite 
conjugare  object  and  image  planes.  Each  object  point  is  linearly  polarized  parallel  to  the 
horizontal  plane.  Since  light'from  each  object  point  is  collimated  before  entering  the 
polarizing  beam  splitter,  every  ray  from  a  given  object  point  encounters  the  same  s-p 
coordinate  orientation.  In  this  situation,  the  polarization  state  from  each  object  point  is 
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transformed  uniformly  by  the  polarizing  beam  splitter.  Figure  4  shows  how  the 
polarization  direction  of  a  set  of  linearly  polarized  object  points  are  rotated  by  the 
polarizing  beam  splitter.  The  intensity  of  each  object  point  is  transmitted  by  the 
polarizing  beam  splitter  according  to 

/(({))  = /oCos2((t))  (2) 

where  (p  is  again  the  rotation  of  the  p  direction  from  the  horizontal  plane,  and  the 
direction  cosines  specify  a  plane  wave  direction.  The  polarizing  beam  splitter  introduces  a 
nonuniform  intensity  distribution  in  the  image  plane  and  variation  of  polarization  direction 
acro.ss  the  image. 

The  transmission  and  reflection  properties  of  a  polarizing  beam  splitter  can  be 
written  in  terms  of  its  Mueller  matrbc  which  is  a  function  of  wavelength  and  angle  of 
incidence.  Appendbc  A  derives  the  Mueller  matrices  in  transmission  and  reflection  of  a 
polarizing  beam  splitter  in  terms  of  the  parameters  generally  calculated  by  thin  film 
analysis  programs,  the  s  and  p  transmittance  and  reflectance,  and  the  retardance  in 
transmi.ssion  and  reflection.  These  Mueller  matrices  include  the  effect  of  the  rotation  of 
s-p  coordinates  with  incident  angle  and  may  be  used  to  facilitate  theoretical  comparisons  of 
e.vperimental  results. 

III.  Figures  of  merit  for  a  polarizing  beam  splitter 

In  this  section,  six  figures  of  merit  which  characterize  the  performance  of  polarizing 
beam  splitters  are  defined.  Measurements  of  these  figures  of  merit  are  reported  in  section 
V. 

The  figures  of  merit  defined  in  this  section  are  measurements  which  determine  how  a 
polarizing  beam  splitter  transmits  and  reflects  two  incident  polarization  states.  The  first 
polarization  state  is  linearly  polarized  parallel  to  the  horizontal  x-z  plane  and  the  second  is 


parallel  to  the  vertical  y-z  plane  (refer  to  Figure  1).  We  will  refer  to  these  directions  as 
horizontal7i  and  vertical  u.  For  an  axial  ray  along  the  z  axis  the  fi  is  aligned  with  p 
direction  and  the  0  direction  aligns  with  the  s  direction. 

The  figures  of  merit  are  measured  as  functions  of  direction  cosines  defined 

in  Figure  1.  The  measurements  are  shown  in  Figures  6  -  9.  The  figures  of  merit, are: 

1.  Transmission  efficiency  T E  -  me  fraction  of  linearly  polarized  light  in  the  /Idirection 
which  is  transmitted  into  the  /rstate. 

2.  Transmission  couplance  TC-  the  fraction  of  incident  linearly  polarized  light  in  the  Fl 
direction  that  is  transmitted  into  theostate, 

3.  Transmission  leakage  TL  -  the  fraction  of  linearly  polarized  light  in  the  o direction 
that  is  transmitted. 

4.  Reflection  efficiency  A?  £-  the  fraction  of  linearly  polarized  light  in  the  c- direction  that 
is  transmitted  into  the  ostate. 

5.  Reflection  couplance  RC -  the  fraction  of  linearly  polarized  light  in  the  ^'direction 
that  is  reflected  into  the  /zstate. 

().  Reflection  leakage  /?/;-  the  fraction  of  linearly  polarized  light  in  the /t  direction  that 
is  reflected. 

Transmi.ssion  and  reflection  efficiency  represents  the  fraction  of  light  that  takes  the 
correct  path  at  the  polarizing  beam  splitter  and  remains  in  the  correct  polarization  state. 
Transmission  and  reflection  leakage  gives  the  fraction  of  light  that  takes  the  incorrect  path 
at  the  polarizing  beam  splitter.  Transmission  and  reflection  couplance  is  the  fraction  of 
light  that  takes  the  correct  path,  but  is  coupled  into  the  orthogonal  polarization  direction. 

The  complete  experimental  polarization  description  of  a  polarization  element  is  in 
terms  of  the  Mueller  Matrix.  The  Mueller  matrix  for  a  polarizing  beam  splitter  in 
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transmission  is  ,vi  < ,  and  the  Mueller  matrix  for  eflection  is  ,  both  of  which  are  functions 
of  wavelength  and  angle  of  incidence. 
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(3) 

Although  the  Mueller  matrices  provide  a  complete  description,  they  are  not  necessarily  the 
clearest  or  most  appropriate  description  for  all  purposes.  For  polarizing  beam  splittersj  the 
six  figures  of  merit  are  readily  measured  and  convey  the  most  important  information.  The 
relationship  of  each  figure  of  merit  to  the  corresponding  Mueller  matrix  elements  is  given 
in  Table  1.  The  relationship  of  the  figure  of  merit  to  the  complex  amplitude  transmission 
and  reflection  coefficients  are  also  given  in  Table  1.  The  corresponding  deri’-ations  are 
given  in  Appendix  A. 

IV.  Imaging  poinrimeter  configuration  and  measurement  procedure 


Figure  3  shows  the  configuration  of  an  imaging  polarimeter  used  to  measure  the 
angular  dependence  of  a  polarizing  beam  splitter.  The  light  source  is  a  850  nm  wavelength 
diode  laser.  The  beam  is  focused  on  a  rotating  ground  glass  disk  to  scramble  the  temporal 
coherence  of  the  beam.  The  beam  is  collimated  and  passes  through  a  silver  glass  polarizer 
(the  polarization  state  generator)  which  is  mounted  in  a  rotary  stage.  A  lens  prepares  a 
converging  spherical  wave  with  a  10°  half  angle  and  is  incident  on  the  polarizing  beam 
splitter.  The  beam  transmitted  (reflected)  by  the  polarizing  beam  splitter  ib  collimated  and 
pa.sses  through  a  second  silver  glass  polarizer  (polarization  state  analyzer).  By  moving  the 
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polarizer/camera  the  same  polarizer  and  camera  are  used  to  make  measurements  both  in 
transmission  and  reflection. 

The  aperture  stop  (see  Figure  3)  of  the  system  is  imaged  onto  the  ccd.  Accordingly, 
each  pixel  represents  an  angular  ray  path  though  the  polarizing  beam  splitter  with  direction 
cosines  (/ ,  m).  By  measuring  the  polarizing  properties  on  a  pixel  by  pixel  basis  the 
polarizing  beam  splitter  is  characterized  as  a  function  of  the  direction  cosines  of  incident 
light. 

Each  image  is  converted  into  a  reduced  resolution  data  set  by  collapsing  4x5  pixel 
regions  into  single  pixels  (super-pixels).  The  unequal  dimensions  of  the  averaged  pixels 
compensate  for  the  5:4  aspect  ratio  of  the  ccd.  Super-pixels  also  reduce  the  intensity 
variation  due  to  speckle  on  the  ccd  and  decrease  the  effect  of  beam  wander  due  to  any 
wedge  (prism)  in  the  rotating  polarization  elements.  High  spatial  resolution  is  not 
important  for  these  measurements  since  the  field  of  view  dependence  of  polarizing  beam 
.splitter  cubes  varies  slowly  with  angle  of  incidence. 

To  control  system  operation,  capture  and  process  images,  and  display  the  results,  a 
PC  386  computer  is  used  with  the  following  peripherals:  a  512  X  512  video  frame  grabber 
board  with  8-bit  resolution  in  conjunction  with  a  visible/near  infrared  monochrome  ccd 
camera.  Data  reduction  is  performed  with  high  precision  floating  point  processing  in  PC 
memory.  A  color  monitor  is  used  to  display  measured  frames  in  false  color. 


V  Measurements  of  the  polarizing  beam  splitter 

Six  figures  of  merit,  described  in  section  III,  were  made  on  four  different  polarizing 
beam  splitters.  The  plots  in  Figures  6-9  are  measures  of  the  angular  dependence  over  a 
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±10°  field  of  view.  Figures  6-8  contain  measurements  made  on  three  different 
commercially  available  polarizing  beam  splitters  from  Melles-Griot,  Newport  and  CVI, 
Measurements  on  a  polarizing  beam  splitter  designed  at  AT&T  Bell  Labs  specifically  to 
operate  over  a  ±7°  field  of  view  is  shown  in  Figure  9. 

Each  figure  of  merit  involves  two  measurements.  In  the  first  measurement,  the 
polarization  analyzer  and  generator  are  oriented  with  their  transmission  axes  along  the 
horizontal  x  or  vertical  y  directions  (refer  to  Figure  1),  according  to  which  figure  of  merit  is 
to  be  measured.  An  image  of  the  exit  pupil  is  digitized  and  stored.  The  image  is  then 
divided  by  the;total  incident  intensity  to  obtain  the  fraction  of  light  transmitted  or  reflected 
from  the  polarizing  beam  splitter  cube. 

The  coordinate  system  used  in  these  plots  is  common  to  both  the  plots  and  the 
polarizing  beam  splitter  shown  in  Figure  1.  Direction  cosine  a  1  increases  from  90®  1 
in  the  positive  x  direction,  f3 1  increases  from  90®  1  in  the  positive  y  direction  and  y  1 
increases  from  0  ®  I  radially.  The  large  circles  indicate  the  zones  where  the  angles  of 
incidence  are  7  =  5°  and  y-  10®.  The  data  is  plotted  as  functions  of  angle  of  incidence  on 
the  entrance  face  of  the  cube  (in  air). 

The  polarizing  beam  splitters  were  measured  with  incident  horizontal  and  vertical 
polarized  light,  i.e,  polarized  in  the  x-z  and  y-z  planes  (refer  to  Figure  1).  Along  the  x  axis, 
tlie,s-p  coordinates  align  with  the  vertical  and  horizontal  directions  i.e.,  along  the  x  axis  the 
s  and  p  transmittances  and  s  and  p  reflectances  (as.t,cip.i,<^s.r,ctp.r)  of  the  multi-layer 
are  represented.  The  positive  x  axis  corresponds  to  an  increase  in  angle  of  incidence  a 
beam  makes  with  the  normal  to  the  beam  splitting  interface,  as  large  as  55®  in  air,  in  the 
positive  X  direction,  and  as  small  as  35®  in  the  negative  x  direction.  We  will  now 
concentrate  our  discussion  along  the  x  axis. 
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The  polarizing  beam  splitters  exhibited  much  wider  angular  bandwidth  for  s 
polarization  than  p  polarization.  The  angular  response  of  p  polarized  light  varied  widely 
among  the  measured  polarizing  beam  splitters.  Here  angles  of  incidence  will  be,given 
relative  to  the  entrance  face  normal,  positive  angles  increasing  in  the  positive  x  direction. 
The  Melles-Griot  cube  exhibited  three  p  state  transmission  peaks  at  -10®,  -2°  and  9°  (Figure 
6a).  Multiple  m^ima  are  indicative  of  field  widening.  Transmission  falls  to  under  80%  at 
5°.  The  s  reflectance  was  very  constant  over  the  measured  angular  range,  about  96% 

(Figure  6b).  A  p  reflectance  maximum  is  found  at  5“  (Figure  6c),  corresponding  to  the 
decrease  in  p  transmission  at  5°.  The  s  transmittance  (Figure  6d)  was  below  1%  for  all 
angles  of  incidence.  As  expected,  coupling  between  s  and  p  polarizations  was  very  small, 
below  1%  (Figures  6e,  6f).  Coupling  between  s  and  p  not  expected  since  the  beamsplitting 
interface  is  isotropic  and  homogeneous. 

The  Newport  polarizing  beam  splitter  exhibited  a  flat  angular  response  to  p 
polarization  for  positive  angles  of  incidence  in  transmission,  but  the  p  transmittance  fell  to 
22%  at  - 10°  (Figure  7a).  The  s  reflectance  of  was  96%  from  -10°  to  10°  (Figure  7b).  A 
sharp  increase  in  p  reflectance  occurs  in  the  negative  direction  (Figure  7c),  explaining  the 
sharp  decrease  off  observed  in  transmission.  S  transmittance  was.small  and  angle 
insensitive  (Figure  7d).  The  couplance  between  s  and  p  polarizations  is  below  1%  (Figure 
7e,  7f). 

The  p  transmittance  of  the  CVI  polarizing  beam  splitter  peaked  at  91%  and 
decreases  rapidly  in  both  directions  (Figure  7a).  The  s  reflectance  was  nearly  98%  over  the 
range  of  angle  of  incidence  (Figure  7b).  The  p  reflectance  increases  correspond  to  the  p 
transmittance  decreases  (Figure  7c).  S  transmittance  was  under  1%  for  all  angles  of 
incidence  (Figure  7d).  Couplance  between  s  and  p  is  below  1%. 
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The  AT&T  polarizing  beam  splitter  was  designed  to.  operate,  over  a  ±1°  field  of  view. 

P  transmittance  was  greater  than  96%  form  -10°  to  T.  The  s  reflectance  varied  from  97% 
at  -10°  to  91%  at  10°.  The  p  reflectance  exhibits  an  increase  from  1%  at  5°  to  7%  at  10°, 
and  is  under  1%  everywhere  else.  The  s  transmittance  increases  from  2%  at  5°  to  6%  at 
10°.  The  couplance  between  s  and  p  is  under  1%. 

The  orientation  of  s  and  p  directions  change  at  the  beam  splitting  interface  for  angles 
of  incidence  in  which  y*  0.  This  property,  explained  in  detail  in  section  II,  is  inherent  to  the 
geometry  of  the  polarizing  beam  splitter  cube  and  degrades  the  performance  in  several 
ways.  First,  the  efficiency  at  which  the  polarizing  beam  splitters  transmitted  horizontally 
polarized  light  and  reflected  vertically  polarized  light  was  decreased  (Figures  6-9,  a,b). 
Second,  the  leakage  of  horizontally  polarized  light  in  reflection  and  vertically  polarized 
light  in  transmission  increased  (Figures  6-9,  c,d).  Finally,  transmitted  horizontally 
polarized  light  was  coupled  into  vertical  and  reflected  vertically  polarized  light  was  coupled 
into  horizontal  (Figures  6-9,  e,f). 

The  polarization  purity  of  transmitted  and  reflected  beams  is  influenced  by  both 
transmittance  and  reflectance  variations  of  the  beamsplitting  interface  and  the  varying 
orientation  of  s-p  coordinates.  For  example,  for  the  Melles-Griot  polarizing  beam  splitter 
the  purity  of  transmitted  horizontally  polarized  light  variations  for  horizontally  polarized 
light  incident,  the  purity  of  transmitted  horizontally  polarized  light  varies  between  160:1  to 
130:1.  The  orientation  variation  causes  a  dramatic  drop  in  polarization  purity  down  to 
30:1.  In  reflection,  the  polarization  purity  along  the  x  axis  remains  constant  at  about  300:1, 
but  drops  to  20:1  due  to  the  rotation  of  s-p  directions. 

The  accuracy  of  our  measurements  are  limited  by  the  A-D  conversion  of  the  8-bit 
camera  used.  VVe  have  not  attempted  to  verify  our  measurements  theoretically  from  the 


14 


equations  in  Appendix  as.the.thin  film.designs  of-the  polarizing  beam  splitters  are  not 
available  to  us.  The  measurements  serve  to  show  the  over-all  angular  behavior  of  the  cube 
and  to  demonstrate  the  spatial  variation  of  linear  polarization  inherent  to  polarizing  beam 
splitters.  The  occasional  abrupt  changes  in  the  data  may  be  due  to  coating  thickness 
variations  and\or  experimental  error. 

V.  Summary 

An  examination  of  several  polarizing  beam  splitter  cubes  has  shown  that  all  are  far 
from  their  ideal  performance.  The  polarizing  beam  splitter  cube  is  intended  to  divide  the 
incident  beam  based  on  its  polarization  state.  In  its  usual  orientation  (beamsplitting 
interface  normal  in  the  plane  parallel  to  our  optical  table),  it  should  transmit  the  horizontal 
component  of  the  incident  light,  and  reflect  the  vertical  component.  However,  a  fraction  of 
the  vertically  polarized  light  is  leaked  in  transmission,  a  fraction  of  horizontally  polarized 
light  is  reflected,  and  some  fraction  of  the  light  which  takes  the  intended  path  has  its  plane 
of  polarization  rotated. 

The  shortcomings  of  any  polarizing  beam  splitter  has  its  source  in  several  places,  (1) 
the  inherent  geometry  of  the  cube,  (2)  the  transmission  and  reflection  coefficients  for  the 
interface,  and  (3)  fabrication  and  tolerancing  difficulties. 

The  inherent  geometry  of  the  polarizing  beam  splitter  cube  arises  because  the 
polarizing  beam  splitter  performs  its  task  in  the  local  s-p  coordinate  system,  and  not  in  the 
global  x-y  coordinate  system.  If  we  were  to  illuminate  it  with  a  beam  with  the  correct 
spatial  variation  of  polarization,  the  ideal  cube  would  transmit  all  of  the  p-component,  and 
reflect  all  the  s-component,  and  divide  our  beams  exactly  as  desired. 
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The  second  shortcoming  arises  because  the  reflecting  and  transmitting  properties  of 
the  interface,  its  Fresnel  coefficients  or  amplitude  transmission  relations,  vary  with  angle  of 
incidence  and  wavelength.  There  can  only  be  a  limited  angular  bandwidth  and  wavelength 
bandwidth  over  which  acceptable  performance  is  achieved.  In  our  data  sets,  the  variations 
of  the  Fresnel  coefficients  with  angle  of  incidence  is  clearly  measured  along  the  x-axis  of 
our  graphs.  As  the  angle  of  incidence  decreases,  the  performance  tends  to  decrease 
rapidly,  as  it  must.  At  zero  angle  of  incidence,  a  homogeneous  and  isotropic  thin  film 
coating  cannot  function  as  a  polarizing  beamsplitting  interface.  Performance  tends  to  stay 
good  for  much  further  toward  increasing  angle  of  incidence.  Thus  many  of  the  polarizing 
beam  splitters  are  improved  if  we  tilt  them  slightly  to  increase  the  internal  angle  of 
incidence. 

We  see  in  the  data  the  great  pains  taken  by  the  thin  film  designers  to  increase  the 
fields  of  view.  Several  cubes  have  Fresnel  coefficients  which  dip,  then  come  back  up, 
indicating  that  they  have  been  optimized  to  have  several  peaks.  This  increases  the  overall 
field,  at  the  expense  of  performance  near  the  center.  (If  the  maxima  were  moved  closer 
together,  it  would  operate  better  over  a  small  field,  than  with  the  maxima  further  apart.) 

Appendix  A:  Mueller  matrix  of  a  polarizing  beam  splitter 

Here  we  derive  the  Mueller  matrbc  for  a  non-depolarizing,  polarizing  beam  splitter 
and  use  the  result  to  find  the  figures  of  merit, listed  in  Table  1  in  terms  of  the  complex 
amplitude  transmission  and  reflection  coefficients  of  the  beamsplitting  interface. 

The  effect  of  an  optical  system  (element)  on  the  polarization  of  incident  light  is  most 
generally  expressed  by 


S  =MS 


(/l.l) 
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where  s  and  s '  are  the  4  x  1  Stokes  vectors  o£the  incident  .and  transmitted  light, 
respectively,  and  iVI  is  the  4  x  4  real  Mueller  matrix  that  characterizes  the  linear  interaction 
of  the  light  with  a  polarizing  optical  element./2<5/  The  elements  of  the  Stokes  vector  are 
defined  as 


A.\ 

Si 

2<£;,><£„>cosy 

\2<£;,><£„>sinY/ 

(/1 .2) 


where  b'  h  and  i'l.are  projections  of  the  electric  field  on  the  /z  and  axes  and  y  is  the 
difference  in  phase  between  the  /zand  t5 polarizations.  The  Mueller  matrbc  of  an  optical 
system  (element)  is  a  function  of  wavelength  and  of  the  angles  of  the  incident  light.  The 
response  of  an  optical  system  to  a  collection  of  spherical  waves  of  arbitrary  polarization  is 
generally  treated  on  a  ray  by  ray  basis  and  can  be  written  as  a  Mueller  matrix  M  ( 2/ ,  P  A) 
as  a  function  of  object  coordinate  H,  pupil  coordinate  p,  and  wavelength  'K.[17,18] 

Two  separate  matrices  are  needed  to  describe  field  of  view  polarization  properties  of 
a  polarizing  beam  splitter  cube:  M < describes  the  beam  splitter  in  transmission,  and  Mr 
characterizes  it  in  reflection.  The  Mueller  matrices  in  transmission  and  reflection  for  an 
ideal  polarizing  beam  splitter  illuminated  at  normal  incidence  are: 
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0 

f  1 

- 1 

0 

1 

1 

0 
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- 1 
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oj 

V  0 

0 

0 
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These  are  the  Mueller  matrices  for  a  horizontal  linear  polarizer  and  a  vertical  linear 
polarizer/id/,  so  the  polarizing  beam  splitter  transmits  horizontally  polarized  light 
S  =  ( 1 .  1.0.0)  and  reflects  vertically  polarized  light  S  =  (1,-1,0,0)l 
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In  transmission,  the  Mueller  matrix  for  a  linear  partial  polarizer,  which  may  also 
display  linear  birefringence  is  [16] 


B 

A 

0 

0 


0 

0 

Ccos6t 
-Csin  6, 


°  1 
Csin6,  r 

CCQSbtJ 


(/1.4) 


where 

5  =  p2  ,(0A)-pL(0A).  (A.S) 

C  ~  2p p_,(0,?v,)Pj_j(0,X,), 

and  p  p . ,  ( 0 ,  \ )  and  p  y . ,  ( 0 ,  X, )  are  the  amplitude  transmittances  for  p  and  s  light,  and  we 
explicitly  consider  that  the  amplitude  transmittance  and  retardance  are  functions  of  angle 
of  incidence  and  wavelength.  6,  =  6p,(  -  6s,„  the  difference  in  accumulated  phase  in 
transmission  between  p  and  s  components,  is  the  linear  retardance  of  the  interface,  0  is  the 
refracted  ray  angle  of  incidence  with  respect  to  the  surface  normal  of  the  beamsplitting 
interface,  given  by 

where  (I .  m .  n)  are  the  incident  ray  direction  cosines  in  air  and  (^ ,  m* ,  n )  are  the 
direction  cosines  in  glass  of  the  ray  after  refracting  at  the  front  face  of  the  beamsplitter 


Q(L,m,n)  =  Arccos 
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cube  given  by 


I  =L, 


m  =m,. 


CA.7) 


n  =n  +  n. 


(/1 .8) 


where  is  the  index  of  refraction  of  the  glass  substrates  and  riois  the  index  of  refraction 
before  the  cube  (generally  air). 

In  reflection,  the  Mueller  matrix  (Eqn.  A4)  is  evaluated  using  the  reflection 


parameters 


/1  =  pJ.,(0.\)-pL(0A). 


fi  =  pf,,(e,x)~p2  ,(e,x), 

^  ~  2Pp,r(0)Pi.r(0*^)* 


(/1.9) 


and  p  p. r  (0  •  ^ )  iind  p  j. (9 ,  A. )  are  the  reflectances  for  p  and  s  light  respectively,  and  6 ,  is 
replaced  with  5^,  the  linear  retardance  between  p  and  s  in  reflection. 

The  incident  and  transmitted  Stokes  vectors  are  defined  in  u  -  h  coordinates.  To 
apply  the  Mueller  matrix  on  the  Stokes  vector,  the  Stokes  vector  must  first  be  rotated  into 
local  s-p  coordinates  of  the  beamsplitting  interface.  This  is  accomplished  by  a  simple 
coordinate  rotation  about  the  normal  to  the  beamsplitting  interface 


where  <|)  is  the  angle  between  horizontal  and  p  directions 

6 ( ^ m , n )  =  Arccosi  -]=•.-  .7  — -(i 

Vv2sin(i>(f,/n,/i)  J 


C/1.10) 


(/1. 11) 


R  (  ({»)  is  the  matrix  for  a  rotational  change  of  basis  for  the  Mueller  calculus/id/. 
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Figure  1  Geometryof  a  polarizing  beam  splitter.  Side  (a,b,c,d)  is  the  entrance  face. 
(a,e,g,d)  is  the  beam  splitting  interface.  The  transmitted  beam  exits  the  face  (f,e,g)  while 
the  reflected  light  exits  at  (b,e,g,c).  An  incident  ray  is  identified  by  its  direction  cosines 
{L,ni,  a)  defined  as  I  =  cos  a,m  =  cos  3 ,  and  n  =  cos  y  •  a  is  the  angle  between  the  x  axis 
and  the  ray.  3  btween  the  y  axis  and  the  ray,  y  between  the  z  axis  and  the  ray.  We  refer  to 
the  x-z  plane  as  the  hprizontal  and  the  y-z  plane  is  vertical.  The  origin  of  the  x-y-z 
coordinate  system  is  located  at  the  center  of  the  entrance  face.  The  cube  is  placed  in  the 
imaging  polarimeter  with  the  horizontal  plane  parallel  to  the  table  top. 


Figure  2  Orientation  of  the  s  and  p  directions  on  the  beamsplitting  interface  of  a 
polarizing  beam  splitter  plotted  as  a  function  of  beam  angle  of  incidence  on  the  cube  face. 
The  orientation  of  the  beam  splitting  interface  is  shown  in  Figure  1.  Figures  1  and  2  share 
the  same  coordinate  system,  ot increases  from  90“  in  the  positive  x  direction.  3  increases 
from  90“  in  the  positive  y,  and  the  large  circles  indicate  the  zones  where  the  angles  of 
incidence  are  y  =  S°  and  y  =  1 0“ .  Incident  rays  for  which  3  >  90°  involve  a 
counterclockwise  rotation  of  s-p  coordinates  and  rays  in  which  3  <  90“  the  rotation  is 
clockwise.  The  polarizing  axis  is  parallel  to  the  p  direction  in  transmission  and  the  s 
direction- in  reflection.  Polarizing  beam  splitters  have  a  spatially  varying  orientation  of 
linear  polarization. 

Figure  3  Polarization  aberration  introduced  into  a  converging  wavefront  focued  through 
an'”ideal  polarizing  beam  splitter.  Incident  on  the  polarizing  beam  splitter  is  a  shperical 
wave  linearly  polarized  in  the  horizontal  x-z  plane.  The  transmission  axis  of  a  polarizing 
beam  splitter  depends  on  ray  angle  of  incidence  and  is  aligned  with  the  p  direction  at  the 
beam  splitting  interface.  The  variation  of  p  direction  with  angle  of  incidence  is  illustrated 
in  Figure  2.  The  incident  polarization  direction  is  shown  in  plot  3a,  and  the  variation  of 
linear  polarization  across  the  transmited  spherical  wave  is  shown  plot  3b.  Variation  of  p 
direction  with  angle  of  incidence  introduces  apodization  and  spatially  varying  linear 
polarization  into  uniformly  polarized  spherical  wavefronts. 

Figure  4  A  ideal  polarizing  beam  splitter  is  situated  in  collimated  space  between 
conjugate  planes.  In  this  situation  a  polarizing  beam  splitter  introduces  intenstity 
non-imiformity  and  polarization  variation  into  the  image.  Light  from  each  object  point  is 
collimated  and  is  incident  at  the  polarizing  beam  splitter  at  different  angle,  the  variation 
of  p  direction  with  angle  of  incidence  (described  in  Figure  2)  introduces  polarization 
direction  variation  in  the  image  plane. 

Figure  5  Configuration  of  the  Imaging  Polarimeter  used  to  measure  four  polarizing  beam 
splitters  at  849nm  over  a  ±  1 0“  field  of  view.  The  source  is  a  849nm  laser  diode.  The 
dashed  line  shows  the  path  of  the  chief  ray.  The  stop  of  the  system  is  imaged  on  to  the 
CCD  detector  and  is  at  the  lens  located  just  after  the  polarizing  beam  splitter  under  test. 
The  polarizing  beam  splitter  is  characterized  as  a  function  of  angle  of  incidence  by 
computing  figures  of  merit  on  a  pixel  by  pwel  basis. 
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Figure  6  Angular  Performance  of  a  Melles-Griot  830  nm  polarizing  beam  splitter 
measured  at  850  nm.  The  plots  on  the  left  are  for  horizontally  polanzed  incident  light,  the 
plots  on  the  right  are  for  vertically  polarized  incident  light. 

(a)  Transmission  efficiency  (R-^  fi)  (b)  Reflection  efficiency  (u  u) 

(cy  Reflection  Leakage  (h /I  +  0  (d)  Transmission  I..eakage  (u  ^  u  + 

(e)  Transmission  Coiiplance  (h u)  (f)  Reflection  Couplance  (u /r) 

(a)  (b)  (c)  (d)  (e)  (f) 

Figure  7  Angular  Performance  of  a  Newport  broad  band  polarizing  beam  splitter 
measured  at  850  nm.  The  plots  on  the  left  are  for  horizontally  polarized  incident  light,  the 
plots  on  the  right  are  for  vertically  polarized  incident  light.  The  plots  are: 

(a)  Transmission  efficiency  (/i (b)  Reflection  efficiency  (u u) 

(c)  Reflection  Leakage  (h-^  R  +  v)  (d)^  Transmission  Leakage  (u  ->  i)  +  /j) 

(e)  Transmission  Couplance  (R  ->  v)  (f)  Reflection  Couplance  {0  R) 

(a)  (b)  (c)-  (d)  (e)  (f) 

Figure  8  Angular  Performance  of  a  CVI  polarizing  beam  splitter  measured  at  850  nm. 
The  plots  on  the  left  are  for  horizontally  polarized  incident  light,  the  plots  on  the  right  are 
for  vertically-polarized  incident  light.  iTie  plots  are: 

(a)  Transmission  efficiency  (R-^  R)  (b)  Reflection  efficiency  {G  0) 

(c)  Reflection  Leakage  {R-^  R  +  v)  (d)  Transmission  Leakage  (v  C  *  R) 

(e)  Transmission  Couplance  (/i  0)  (f)  Reflection  Couplance  (u  R) 

(a)  (b)  (c)  (d)  (e)  (f) 


Figure  9  Angular  Performance  of  a  polarizing  beam  splitter  desinged  at  AT&T  Bell  Labs 
to  operate  over  7°  field  of  view.  Measurements  made  at  850  nm.  The  plots  on  the  left  are 
for  horizontally  polarized  incident  light,  the  plots  on  the  right  are  for  vertically  polarized 
incident  light.  The  plots  are: 

(a)  Transmission  efficiency  (R-^  R)  (b)  Reflection  efficiency  (0  0) 

(c)  Reflection  Leakage  (/t h  +  u)  (d)  Transmission  Leakage  (u u  + 

(e)  Transmission  Couplance  (R  ->  0)  (f)  Reflection  Couplance  (0  R) 

(a)  (b)  (c)  (d)  (e)  (f) 


Figure  of  Merit 


Corresponding  Mueller  matrix  Corresponding  amplitude  coefficient 
expression  expression 


Transmission  efficienev 
rF(a.fJ) 

|(m',,  +  m',2+m^,  +  /n22) 

Pp  ,cos*p-  p;  jSin^p  -  2p  p_,p,_,sin-pcos'Pcos6, 

Transmission  leakage 
77.(a.(J) 

-(/n,|  -  /n  ,2) 

^(pf,.,cos^fi  +  pf  ,Gin^P) 

Transmission  couplance 
VCCa.l^) 

sin^pcos^Pfp^,  ,*  p^  ,-2p„  ,p,  ,cos6,) 

Reflection  efficiency 

m'tz-  mil* 

p*,,cos‘P'  p*_,sin‘p-  2p,,,pp_,siii*|Jcos'pcos6, 

Reflection  Lcakaue 

1  ,  ,  - 

-(p‘  pC0S‘P  +  p;  ,ein-p) 

Rellection  coufilaiice 
A'C(a.li) 

sin^Pcos^Pfp^  P  +  pj;  ,.-2p,  ^p,,  ,.cos6p) 

Table  J;  Measured  Mueller  matrbe  elements  and  theoretical  values  of  angular  performance  parameters 
of  polarizing  beam  splitters.  p(a(0.<|)),\)is  the  real  amplitude  transmission  or  rellection  coefficient. 
Subscripts  p  and  s  refer  to  p  polarization  and  s  polarization  respectively.  Superscripts  t  and  r  denote 
transmi.ssion  and  rellection.  6, and S^are  the  retardance  in  transmission  and  rellection  respectively. 

<r (/./)! .  /i)  is  the  angle  of  incidence  a  ray  makes  with  the  bcamsplitting  interface  and  P(/ .  rn .  /i)is  the 
orientation  of  s-p  coordinates  relative  to  ft  -  v  coordinates.  The  functional  dependence  of  a  ( / .  m ,  n ) 
andfiff  given  in  appendix  A. 


Figure  3 


Figure  4 


Marginal  ray 


Transmission  efficiency  (h  /i) 


Reflection  efficiency  (0  ->  v) 


Transmission  couplance  (h  0)  Reflection  couplance  (v  fi) 


